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I.  INTRODUCTION 


Modern  finite  difference  procedures  for  solving  the  partial  differential 
equations  which  describe  fluid  flow  frequently  utilize  curvilinear  mapping 
procedures.  Because  boundary  surfaces  in  the  physical  plane  can  be  mapped 
onto  rectangular  surfaces  in  the  transformed  plane,  a  finite  difference 
algorithm  for  the  transformed  equations  can  be  readily  applied  to  a  variety  of 
different  body  shapes.  Even  unsteady  body  motion  is  easily  incorporated  into 
the  governing  equations.  To  take  advantage  of  the  generality  of  the  trans¬ 
formed  equations,  however,  one  needs  a  fairly  automatic  method  of  generating 
smoothly  varying  grids  that  fit  arbitrary  bodies  and  allow  grid  point 
clustering.  The  problem  of  grid  generation,  restricted  to  arbitrary  pro¬ 
jectile  shapes,  is  the  subject  of  this  report. 

A  modular  general  purpose  grid  generation  routine  has  been  written  for 
use  with  standard  and  nonconventional  projectiles  shapes.  Three-dimensional 
grid  generation  capability  is  satisfied  by  generating  a  sequence  of  planar 
grids  about  axis  normal  cross  sections.  The  grid  generation  routine  can  be 
used  to  generate  planar  Cartesian-like  grids,  C-grids,  0-grids,  or  any  portion 
thereof.  The  routine  can  be  used  for  axisymmetric  projectiles  with  or  without 
stings,  symmetric  tubular  projectiles,  and  for  any  projectile  with  an  axisym¬ 
metric  nose. 

The  modular  grid  generation  program  developed  here  contains  its  own  body 
surface  representation  and  grid  point  distribution  routines.  Another  set  of 
routines  allows  the  user  to  build  up  an  arbitrary  outer  boundary  curve  and 
grid  point  distribution.  Finally  the  mesh  itself  is  formed  using  either 
algebraic  straight-line  rays  to  connect  inner  and  outer  boundary  points,  or 
by  using  an  elliptic  solver.  Clustering  along  grid  lines  is  accomplished 
with  an  exponential  clustering  routine  that  allows  the  user  to  specify  a  given 
grid  point  spacing  along  the  inner  boundary  curve.  The  modular  structure  of 
the  grid  generation  program  allows  the  user  to  substitute  alternate  clustering 
and  grid  generation  routines  of  his  own  design. 

Figures  1-8  show  the  various  classification  of  grids  which  the  code  is 
designed  to  handle.  These  include  axisymmetric  projectiles  with  or  without  a 
sharp  leading  edge  and  with  or  without  an  afterbody  sting.  Figures  1-3. 
Isolated  boattails  or  flares  can  also  be  meshed  (Figure  2b).  The  code  can 
also  treat  tubular  projectiles  (ring  airfoils)  that  are  axisymmetric.  Figure 
4.  The  isolated  axisymmetric  blunt  body  problem,  Figure  5,  is  actually  a 
special  case  of  Figure  1. 

In  Section  II  the  mechanics  of  the  grid  generation  routine  are  described. 
Various  grids  are  displayed  to  illustrate  the  ideas.  Additional  discussion  of 
the  computer  grids  is  given  in  Section  III,  while  complete  documentation  of 
the  codes  is  given  in  the  Appendices.  The  option  of  selecting  a  2  or  3  dimen¬ 
sional  grid  is  included  and  a  plotting  routine  is  presented  in  the  appendix. 


7 


II.  GRID  GENERATION 


The  purpose  of  the  grid  generation  routine  is  to  generate  a  network  of 

constant  lines  of  £  and  n  in  the  physical  x-y  plane  as  indicated  in  Figure 

9a.  Correspondi ng  uniform  values  of  £  and  n  in  the  computational  space  define 

a  one  to  one  mapping  between  points  j,k  in  the  physical  plane  to  points  j,k  in 

the  computational  plane,  see  Figure  9b.  The  mapping  functions  are  described, 

at  least  numerically,  once  £.  ,  and  n.  ,  are  known  in  the  physical  plane  as  a 

J  j  ,k 

function  of  Xj^  and  yj^;  or  conversely,  once  Xj  ^  and  yj  ^  are  determined  in 

the  transform  plane.  The  metric  quantities  £  ,  £  ,  n  ,  and  n,  needed  in  the 

x  y  x  y 

transformed  flow  equations  can  then  be  determined  numerically  (see,  for 
example.  References  1-3). 

In  Figures  10a  and  10b  we  show  a  typical  grid  generated  for  an  axisym- 
metric  projectile  shape.  In  this  case  a  spherical  cap  is  placed  at  the  end  of 
the  boattail  (see  Figure  10b)  in  order  to  avoid  a  highly  discontinuous  corner. 
The  grid  as  it  stands  is  suitable  for  axisymmetric  or  n-invariant1 2 3 4  flow 
calculations.  A  three-dimensional  grid  can  be  formed  by  rotating  the  grid 
about  the  axis  and  defining  constant  lines  in  0.  When  viewed  in  this  manner 
the  grid  is  equivalent  to  warping  a  spherical  coordinate  about  a  nonspherical 
projectile  shape.  Using  the  notation  defined  in  Figure  9a,  the  axis  corres¬ 
ponds  to  the  £  =  0  and  £  =  £  lines.  The  inner  and  outer  boundaries  corres- 

md  x 

pond  to  n  =  0  and  n  =  n 

max 

Depending  on  the  projectile  shape  a  warped  hemispherical  coordinate  may 
be  used  (Figure  1)  or  a  warped  Cartesian  coordinate  (Figure  3).  A  ring  air¬ 
foil  (tubular  projectile)  can  be  meshed  with  a  C-grid  (Figure  4).  One  could 
also  spin  an  airfoil  0-grid  (Figure  7)  about  the  tubular  projectiles  axis  of 
symmetry.  Whichever  class  of  grid  is  used,  however,  it  must  map  onto  the 
uniform  computational  plane  shown  in  Figure  9b.  The  r,  =  0  plane  is  not  neces¬ 
sarily  restricted  to  only  the  body  surface.  It  may  for  example,  include  the 
forward  cut  of  Figure  2a  or  the  lower  and  upper  cuts  of  Figure  4. 


1.  Steger,  J .  L.,  "Implicit  Finite-Difference  Simulation  of  Flow  About 
Arbitrary  Two-Dimensional  Geometries",  AIAA  Journal,  Vol.  16,  July  1978, 
pp.  679-686. 

2.  Pulliam,  T.  H.  and  Steger,  J.  L.  "On  Implicit  Finite-Difference 
Simulations  of  Three-Dimensional  Flow",  AIAA  Paper  No.  78-10,  1978. 

3.  Schiff,  L.  3.  and  Steger,  J.  L.,  "Numerical  Simulation  of  Steady  Super 
sonic  Viscous  Flow",  AIAA  Paper  No.  79-130,  1979. 

4.  Nietubicz,  C.  J.,  Pulliam,  T.  H.,  and  Steger,  J.  L.,  "Numerical  Solution 
of  the  Azimuthal- Invariant  Thin-Layer  Navier-Stokes  Equations",  AIAA  Paper 
No.  79-0010,  1979. 
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In  generating  a  projectile  grid  such  as  those  indicated  in  Figures  1-8, 
one  first  decides  what  class  of  grid  fits  the  given  problem.  The  grid  genera¬ 
tion  problem  can  then  be  broken  into  three  main  tasks  as  follows:  (1)  define 
the  body  shape,  possible  sting  or  cut,  and  distribute  grid  points  along  the 
n  =  0  boundary  (i.e.,  specify  K  as  a  function  of  x  and  z  along  n  =  0).  Points 
along  this  boundary  should  be  clustered  to  flow  field  gradients,  e.g.,  the 
forward  stagnation  point,  expansions,  shocks;  (2)  define  the  outer  boundary 
curve  and  distribute  grid  points  along  the  n  =  boundary.  Here  we 

restrict  K  =  0  and  £  =  £max  to  be  vertical  or  horizontal  straightline  rays  in 

order  to  simplify  programming  logic,  thus  the  endpoints  of  the  n  .  curve  must 

ITlaX 

properly  align  with  those  of  the  n  =  0  boundary;  (3)  once  the  outer  boundaries 
are  defined,  they  are  "connected"  by  generating  the  interior  grid  with  appro¬ 
priate  clustering  functions  in  n. 

In  the  remainder  of  this  Section  the  procedures  used  to  generate  the 
n  =  0  boundary,  the  n  =  n  boundary,  and  the  interior  clustered  grid,  will 
be  described.  max 

A.  Surface  Representation  and  Grid  Point  Distribution 

The  first  step  in  generating  the  grid  is  to  represent  and  distribute 
points  along  the  body  surface.  A  sting  or  cut  may  also  have  to  be  included. 

The  body  shape  is  expected  to  have  either  an  analytic  description  or  be 
described  as  a  table  of  x,y  ordinates.  In  either  case  the  data  is  assumed  to 
be  nondimensional  with  respect  to  a  reference  length  which  can,  of  course,  be 
taken  as  1  so  the  data  remains  unaltered. 

The  present  code  allows  for  either  a  parabolic  arc  or  standard  class  pro¬ 
jectile,  such  as  a  projectile  with  a  sharp  tangent-ogive  or  blunt  secant  ogive 
nose,  cylindrical  body,  boattail,  and  spherical  cap.  If  the  analytic  body 
shape  differs  from  the  above  mentioned  shapes,  then  the  user  must  supply  his 
own  description.  In  this  case  values  of  x  along  the  body  axis  (or  chord)  will 
be  distributed  by  contiguously  combining  segments  of  the  clustering  function 

xi  =  X  +  a<K  +  b\Jj?  +  Ci|>?  o  j  f  (1) 

J  J  J  J  jQ  <  j  <  jf 

where  =  ( j- jQ)/( jf-jQ)  and  j  is  an  index  value  such  that  points  j0  to  jf 

lie  in  the  interval  x0  to  Xf  and  xj^  =  x0  while  Xj  =  Xf.  Equation  (1)  is 

used  to  cluster  Xj  as  a  function  of  j  as  indicated  in  Figure  11.  The  user 

determines  the  shape  of  the  clustering  function  by  specifying  the  initial  and 
final  increments  of  x,  that  is 
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(2a) 


Ax°  ■  V  -  xj„ 


Vxf  =  xJf  ‘  xJf-l 


Since  x0  and  Xf  are  also  specified,  a,  b,  and  c  are  determined 


(2b) 


c  =  {VXf  +  axq  -  2h(x^-xQ) }/(h  -  3h2  +  2h3) 

b  =  {axq  -  h(xf-xQ)  -  c(h3  -  h) } / ( h 2  -  h) 

a  =  x*  -  x  -b-c 
f  o 

where  h  =  (jf  -  j0)_1. 


The  amount  of  clustering  at  each  point  is  determine  by  the  specified 

values  of  AxQ  and  Vx^.  Moreover,  becajse  axq  and  Vx^  are  specified,  the  user 

can  smoothly  patch  functions  together  to  form  a  general  clustering  function. 
Examples  of  this  are  indicated  in  the  computer  output  presented  in  Appendix 
C.  One  drawback  to  the  clustering  function,  Eq.  (1),  is  that  the  function  is 
not  guaranteed  to  be  monotone  in  the  interval.  This  can  happen,  for  example, 
if  Axq  is  too  small  and  Vx^  too  large.  Again,  the  output  in  Appendix  C  indi¬ 
cates  practical  values  to  choose. 

In  the  case  of  a  nonanalytic  body  shape  x  and  y  are  read  in  as  a  table  of 
values.  Here  either  an  axis  length  or  surface  arc  length  is  used  as  a  clus¬ 
tering  function.  If  the  axis  length  is  chosen,  x  is  given  and  corresponding  y 

values  are  found  from  the  table  of  x,y  coordinates  using  cubic  spline 

interpolation5.  Alternately,  the  surface  arc  length  can  be  used.  In  this 
approach  the  arc  length  s  is  computed  from  the  table  of  x,y  values,  and  the 
length  is  normalized.  A  new  normalized  clustered  arc  length  is  then  defined, 
using  Eq.  (1)  with  s  (the  arc  length)  in  place  of  x.  Both  x  and  y  are  then 
interpolated  from  the  tables  x  versus  s  and  y  versus  s.  Again  cubic  spline 
interpolation  is  used. 

Finally,  a  sting  or  cut  may  be  added  to  the  configuration.  Again,  the 
clustering  relation,  Eq.  (1),  is  used  to  distribute  ooints  along  the  sting  or 
cut. 


5.  Dahlquist,  G.  and  Bjorck,  A.,  "Numerical  Methods".  Prentice-Hall , 
Englewood  Cliffs ,  Hew  Jersey ,  1974. 
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B.  Outer  Boundary  Formation  and  Grid  Point  Distribution 

Usually  the  shape  of  the  outer  boundary  curve  is  not  as  well  defined  as 
the  inner  boundary.  However,  the  shape  of  the  outer  curve  may  be  partially 
prescribed.  For  example,  in  Figure  4,  the  axis  of  symmetry  of  the  ring 
airfoil  has  a  specified  location.  The  side  and  upper  boundaries  need  only  be 
smooth  curves  far- removed  from  the  body.  If,  however,  the  ring  airfoil  or 
standard  projectile  is  tested  in  a  wind  tunnel,  a  numerical  simulation  of  the 
experiment  requires  a  fixed  wall  outer  boundary.  In  this  case  the  top  portion 
of  the  outer  boundary  curve  must  be  specified. 

A  part  of  the  grid  generation  problem  then  is  the  formation  of  the  arbi¬ 
trary  outer  boundary.  Here  this  boundary  is  built  up  by  connecting  contiguous 
cubic  segments,  which  in  the  degenerate  case  can  be  straight  lines.  Figures 
12a  and  12b  illustrate  two  typical  outer  boundary  curves.  In  Figure  12a  three 
cubic  segments  make  up  the  boundary  n  =  nmax*  Each  segment,  from  a  to  b  for 

example,  is  formed  by  specifying  the  endpoints  x,y,  and  angle  0,  where  0  is 
the  angle  between  the  curve  and  the  x  axis.  In  the  example,  Figure  12a, 

0  =  90°,  0  =  0  =  0°  or  180°  and  0.  =  90°. 

a  b  c  d 

The  data  x,y,0  at  each  endpoint  determines  the  shape  of  the  parametric 
curves 


x  =  x  +  a,  t  +  cut2 
o  1  2 


o 


0  <  t  <  1 


(4) 


y  =  yo  +  V  +  e2t2 


which  are  equivalent  to  a  cubic 


y  =  y0  +  y1(x-xq)  +  y2(x-xq)2  +  y3(x-xq)3 


(5) 


dy 

The  parametric  cubic  is  used  because  the  condition  =  °°  can  be  specified, 
e.g.,  segment  be  of  Figure  12b  has  this  constraint  at  both  endpoints. 


The  solution  of  the  parameters  ,  B^ ,  and  B2  are  given  by 


(6) 


e2  "  <*f  -  *o>  ‘  61 


11 


where 


and 


mo  = 


dx 

"31 


n  -  dy 

n  —  —  t  x 
o  at 


o 

o 


dy 

n  =  m 
o  o  dx 

dy 

nl  =  ml  lx 


m 


dx 

1  ~3t 


nl 


1 

1 


or 


_  _  _  dx 

m  _  n  — — 
o  o  ay 


or 


_  _  n  dX 
m,  =  n,  — r- 

1  1  dy 


The  solutions  for  nQ,  mQ,  nj,  and  mj  are  conditional  insofar  that  infinite- 
slopes  are  avoided.  Then  the  regular  solutions  are 


i)  if  yxl0  ^  xylo»  yxli  ^  xy 1 1 »  an<i  yx 1 1  *  ^xlo 

mi  =  2  [(Xf  -  y0)  -  (xf  -  x0)  yx|03/(yxli  -  yxl0) 
m0  -  2  ( xf  -  x0)  - 

(7a) 

no  mo  yxio 

ni  =  mi  yxli 

ii)  If  yx|0  <  xy|0,  xy|j  <  yx|j,  and  xy | Xyx | 0  *  1 

mo  -  2^-(xf  ■  x0)  _  (yf  “  y0)  xy I i 3 / ( i  -  xy I i  yxl0) 

nl  =  2  (yf  -  y0)  -  m0  yx|0 

(7b) 

no  -  mo  yx  I  o 
mi  =  nj  xy|j 

iii)  If  xy | q  <  yx|0>  yx 1 1  <  xyli»  an(I  y x 1 1  xy I o  *  I 

no  "  2C(yf  -  y q)  -  (xf  -  xQ)  yx I iH/(i  -  yxli  xy|0) 

ml  “  2  (xf  “  xo^  -  n0  xy I o 

(7c) 

mo  "  no  xy I o 

nl  =  ml  yxli 
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(7  d) 


1111)  If  xy | o  <  yxl0’  xy 1 1  <  yxll»  and  xy 1 1  *  xy I o 

Hi  =  2  [(xf  -  Xq)  -  (y-f  -  y q)  Xy  I  o*) / ( Xy  1 1  “  xy  I  o) 

no  “  2  (yf  -  y0)  -  nl 

m0  Hq  Xy | q 

ml  =  nl  xy  1 1 


Whenever  the  third  constraint  is  violated  (for  example  in  case  (i),  if 
yx | i  =  yx|0)  3  linear  curve  is  used.  In  this  case  ^  =  0  and 

aj  =  x^-x  ,  3^  =  yf-y0.  The  segments  ab,  be,  and  ea  of  Figure  12b  are  exam¬ 
ples  of  the  straight  line  segments. 

The  outer  boundary  curve  is  made  up  of  contiguous  cubic  segments  starting 
from  the  5=0  boundary.  Points  are  distributed  along  this  curve  either  as  a 
uniform  distribution  of  arc  length,  or  as  a  specified  arc  length  distribution 
using  the  previously  defined  clustering  scheme,  Eq.  (1).  Since  the  true  arc 
length  is  not  specified  a  priori,  precise  alignment  of  points  along  the  outer 
boundary  can  be  specified  only  after  the  cubic  segments  are  specified  and  the 
arc  length  is  computed.  Without  such  knowledge  a  normalized  clustering 
function  should  be  used. 


C.  Grid  Generation  and  Clustering 

The  task  of  generating  a  grid  is  undertaken  once  the  boundary  curves  are 
specified  and  points  are  distributed  on  the  n  =  0  and  nmax  boundaries.  Two 

types  of  grid  generation  procedures  were  used  and  are  discussed  below. 

In  the  first  case,  lines  of  constant  5  (i.e.,  the  rays  emerging  from  the 
body)  are  formed  by  simply  connecting  straight  lines  from  points  along  n  =  0 
to  points  along  n  =  nmax«  The  spacing  in  n  along  each  such  line  is  either 

uniform  or  is  determined  by  the  relation 


Ask  -  AS0(1  ♦  e)k->  ,  k  .  1,  kmax  -1  (8) 

Here  Asq  is  the  specified  constant  grid  spacing  at  the  inner  boundary.  The 

parameter  e  is  determined  by  a  Newton-Raphson  iteration  process  so  that  the 

sum  of  the  above  increments  matches  the  known  arc  length  between  the 

n  =  0  and  n  =  n  ,  for  points  which  have  the  same  values  of  5.  Figures  13a 
max 

and  13b  illustrate  a  straight  ray  grid  with  clustering  in  n  for  a  tubular 
projecti le. 
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In  the  second  case,  the  grid  is  generated  with  elliptic  partial  differ¬ 
ential  equations  following  References  6,  7,  8.  The  grid  generating  equations 
are  solved  on  the  specified  computational  space  for  unknowns  Xj  ^  and  y,-  ^ : 


where 


ax  -  2bxc  +  yx  =  -o  (Px.  +  Qx  ) 
5S  ?n  nn  v  t  n' 

ay__  -  2By,.  +  yy  =  -0  (Py_  +  Qy  ) 

JZE,  J?n  Jnn  v  Jn; 


(9a) 

(9b) 


2  2 
x  +  y 
n  n 


3  = 


XpX  +  y  y 
5  n  J  K  n ' 


2  A  2 
=  XC  +  ’ 


J  = 


XrY  - 

KJ  n 


X  Yr 

rr  5 


P  =  P  e~d(r|"n0)  +  P  e"a^n"r)max^ 
o  m 

7T  n  -b(n-n  )  ,  r,  -b(n-n  ) 

Q  =  Q„  e  '  o'  +  Q  e  '  max' 
o  m 


Here  PQ,  Q0,  Pm,  Qm,  a  and  b  are  prescribed  clustering  parameters.  Along 
the  n  =  0  and  n  =  nmax  boundaries,  Xj^  and  yj^  have  been  previously 
prescribed.  Along  the  £  =  0  and  £  =  £  ,  which  are  either  vertical  or  hori- 

Hla  X 

zontal  lines  in  the  physical  space,  the  following  boundary  conditions  are 
enforced:  either 


x  is  given  and  y^  =  0 

(10a) 

on  a  vertical  boundary,  or 

x^.  =  0  and  y  is  given 

(10b) 

on  a  horizontal  boundary. 

6.  Chu ,  V.  H.,  "Development  of  a  General  Finite  Difference  Approximation  for 
a  General  Domain".  Journal  of  Comp.  Physios ,  Vol.  8,  1971 ,  pp.  392-408. 

7.  Thompson ,  J.  F.,  Thames ,  F.  C.,  and  Mastin,  C.  M.,  "Automatic  Numerical 
Generation  of  Body-Fitted  Curvilinear  Coordinate  System  for  Field 
Containing  any  Number  of  Arbitrary  Two-Dimensional  Bodies".  Journal  of 
Comp.  Physics ,  Vol.  15,  1974 ,  pp.  299-319. 

8.  Sorenson,  R.  L.  and  Steger,  J.  L.,  ", Simplified  Clustering  of  Nonorthogonal 
Grids  Generated  by  Elliptic  Partial  Differential  Equations" .  NASA 

TM  73252,  August  1977. 
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For  periodic  grids  as  indicated  in  Figure  7,  these  boundary  conditions  in  £ 
are  replaced  by  the  usual  periodic  relations. 

The  derivative  expressions  on  the  left  hand  side  of  Eq.  (9)  are  all 
differenced  with  conventional  second  order  central  difference  operators,  that 
i  s 


x 

x 


5 

n 


x 

x 

X 


nn 


<xj+i,k  - 
<xj,k+l  ' 
(xj+l,k  - 
(xj+l,k+l 
(xj,k+l  ' 


Xj.l>k)/(2A5) 

xj,k-l>/<21"> 


xj-l,k>/<45>2 


"  Xj+l,k-l 


"  XJ-l,k+l 


+ 


xj-l  }k_i)/(4A^An) 


(11) 


while  derivatives  of  y  are  treated  identically.  The  Jacobian  J  is  formed  with 

central  differencing.  The  right  hand  side  companion  terms  to  P  and  Q, 

_howeve_r,  are  _back ward  or  forward  differenced  depending  on  the  sign  of 
P  and  Q.  If  P  is  positive,  and  y^  are  forward  differenced.  The  terms 

x  ,  y  are  differenced  in  the  same  way. 
n  wn  J 

The  one  sided  differencing  for  the  right  side  term  was  chosen  assuming  J 
is  a  constant.  Preliminary  analysis  with  local  linearization  of  terms  like 
J2x^  suggests  one  sided  differencing  should  also  be  used  in  J  to  keep  balanced 

coefficients.  This  however  has  not  been  evaluated. 


The  difference  equations  to  Eq.  (9)  are  solved  with  a  successive  line 
overrelaxation  (SLOR)  procedure.  As  an  initial  guess  for  the  relaxation 
procedure  we  use  the  straight  Jjne  rjy  procedure  previously  described.  For 
the  most  part,  if  coefficients  P  and  Q  are  large,  the  SLOR  procedure  is  very 
difficult  to  converge.  Consequently,  we  recommend  using  the  algebraic  clus¬ 
tering  function,  Eq.  (8). 

In  the  algebraic  clustering  approach  the  elliptic  solver  is  used  to 
generate  a  grid  with  P  =  Q  =  0.  The  x,y  points  along  a  5  =  constant  line  are 
then  redistributed  along  this  line  as  a  function  of  arc  length.  The  cluster¬ 
ing  function  Eq.  (8)  is  used  for  this  purpose.  This  procedure  works  quite 
well  and  provides  excellent  control  of  the  grid  spacing  near  the  body  surface. 
Further  details  are  given  in  Reference  8.  The  grid  shown  in  Figure  14  was 
generated  in  this  manner. 

The  elliptic  solver  need  not  be  used  over  the  entire  range  in  5.  Because 
of  the  boundary  condition,  Eq.  (10),  the  elliptic  equations  can  be  joined  to  a 
straight  ray  along  any  vertical  or  horizontal  boundary  line  in  5.  Figure  15 
shows  details  of  such  a  procedure  used  in  the  previous  tubular  projectile 
case.  Here  the  5-region  over  the  tubular  projectile  is  meshed  using  the 
elliptic  equations  while  the  remainder  is  meshed  with  straight  rays.  After 
the  basic  grid  is  formed,  the  entire  grid  is  clustered  in  n  using  Eq.  (8). 


15 


D.  Grid  Plotting 


An  integral  part  of  the  grid  development  program  is  the  ability  to  plot 
the  computed  grid  in  a  timely  manner.  The  plot  program  which  was  developed 
and  utilized  allowed  almost  instantaneous  viewing  of  the  computed  grid.  This 
capability  significantly  reduced  the  grid  generator  development  time. 

The  plot  program  was  written  using  Tektronix  Plot  10  software  on  the  BRL 
Cyber  176  computer.  A  program  listing  is  presented  in  Appendix  D. 

The  only  input  required  for  the  plotting  program  is  the  converged  grid 
file  and  the  minimum  and  maximum  x,y  values  of  the  grid.  The  interactive 
program  uses  prompts  for  the  remaining  input. 


III.  DISCUSSION  OF  RESULTS  AND  CONCLUDING  REMARKS 

Figures  10,  13,  14,  and  15  give  the  reader  a  reasonably  clear  picture  of 
the  capability  of  the  grid  generation  routine.  The  other  grids  classified  in 
Figures  1-8  simply  use  the  same  program  elements  in  different  arrangements. 

An  analytic  shape  that  was  meshed  using  the  elliptic  equation  approach  is 
illustrated  in  Figure  10.  Along  the  body,  points  are  clustered  to  the  nose, 
boattail  junctures,  and  base.  No  cuts  or  stings  are  used.  Only  two  cubic 
segments  are  used  to  define  the  outer  boundary  and  along  this  curve  points  are 
uniformly  distributed.  Solution  of  Eq.  (9)  with  no  additional  reclustering 
completes  the  grid  generation  problem.  Note  that  Eq.  (10b)  is  well  satisfied 

along  the  5=0  and  5=5  axis. 

3  max 

In  Figure  13  a  grid  for  a  tubular  projectile  is  shown.  The  body  is 
defined  by  x,y  ordinates  and  here  upper  and  lower  cuts  are  used.  The  outer 
boundary  is  defined  using  four  cubic  segments,  two  of  which  degenerate  to 
straight  lines.  From  the  trailing  edge  on  back,  the  point  distribution  along 
the  outer  boundary  matches  that  of  the  cuts.  In  this  way  vertical  rays  are 
used  over  the  cut,  although  this  is  not  required.  Straight  line  rays  make  up 
the  interior  grid,  and  along  these  rays  points  in  n  are  exponentially  cluster¬ 
ed  using  Eq.  (8).  The  controlled  grid  spacing  along  the  body  is  illustrated 
in  Figure  13b. 

The  case  shown  in  Figure  14  is  similar  to  Figure  10  only  here  the  grid 
generated  using  Eq.  (9)  was  reclustered  along  lines  of  constant  5.  The  grid 
spacing  near  the  body  is  now  controlled  as  before. 

Finally,  the  case  shown  in  Figure  15  is  similar  to  that  of  Figure  13  only 
now  an  elliptic  solver  is  used  over  the  airfoil.  The  cut  region  is  again 
treated  with  vertical  rays. 

Grids  for  nonaxi symmetric  bodies  with  axisymmetric  noses.  Figure  6,  can 
be  generated  as  follows.  For  the  axisymmetric  nose.  Figure  5,  the  grid  is 
generated  in  a  plane  and  then  spun  around  the  axis  forming  a  three-dimensional 
grid.  The  remaining  grid  can  then  be  generated  by  taking  planar  cuts  normal 
to  the  axis  at  various  increments  Ax  (x  aligned  with  the  axis,  see  Figure 
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6).  At  each  cut  a  planar  grid  is  generated  and  the  combinations  of  these 
grids  form  the  three-dimensional  mesh.  At  each  cross  section  one  generates 
the  0-type  grids  shown  in  Figures  7  and  8  being  careful  to  maintain  continuity 
in  x. 


Completed  computer  code  documentation  is  provided  in  Appendices  A  and 
B.  Input  and  output  to  obtain  the  grid  shown  in  Figure  15  is  included  in 
Appendix  C  and  the  plotting  code  is  given  in  Appendix  D. 

The  modular  program  developed  here  has  proven  to  be  quite  flexible,  and 
should  find  application  in  determining  grids  for  various  conventional  and 
nonconventional  projectile  shapes. 
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Figure  1.  Standard  Projectile  Grid  with  Sting 


Figure  2a.  Cartesian-like  Projectile  Grid  with  Sting 
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Figure  2b.  Special  Case  Isolated  Boattail 
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Figure  3.  Standard  Projectile  Grid  with  Base, 
Special  Case  of  0-Grid  with  Symmetry 
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Figure  4.  Tubular  Projectile  Grid  or  C-Grid 
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Figure  5.  Projectile  Blunt  Body  Grid  (Fraction  of  0-Grid) 


Figure  6.  Projectile  with  Symmetric  Nose  and  Nonsymmetric  Afterbody 
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Figure  9.  Mapping  from  Physical  Space  to  Computational  Space 
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Figure  10a.  Standard  Projectile  Grid  in  the  Physical  Plane 
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Figure  10b.  Grid  Detail  near  Projectile 


Figure  11.  Stretching  Function,  Points  j  are  Specified  along 
with  xQ,  xf,  axq,  and  vxf 
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(a)  STANDARD  PROJECTILE  GRID 
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Figure  12.  Outer  Boundary  Structure  and  Terminology 
for  Two  Classes  of  Grid 
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Figure  13a.  Overview  of  Straight  Ray  Grid  for  Tubular  Projectile 
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Figure  13b.  Grid  Detail  near  Lower  Trailing  Edge 
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Figure  14a.  Standard  Projectile  Grid  with  Controlled  Reclustering  along  Lines  of  Constant  £ 
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Figure  14b.  Boattail  Base  Detail  Showing  Specified  Clustering  in  n  near  the  Body 
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Figure  15.  Hybrid  Straight  Ray  and  Elliptic  Solver  Grid  Detail  for  a  Tubular  Projectile 
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Jacobian  of  the  transformation  between  the  physical  and  the  compu¬ 
tational  coordinates 

clustering  parameters  for  the  elliptical  solver 
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APPENDIX  A 


COMPUTER  CODE  DESCRIPTION  AND  LISTING 


The  computer  program  is  a  highly  modular  code.  The  main  program  is 
divided  into  three  parts:  inner  boundary,  outer  boundary,  and  grid  genera¬ 
tion. 

A.  Inner  Boundary 

In  forming  the  inner  boundary,  subroutine  BODY  is  called  to  define  the 
body  shape  and  to  distribute  points  along  the  body.  Subroutine  BODAN  is 
called  for  an  analytic  shape.  The  user  can  modify  this  routine  to  supply  his 
own  body  function.  Subroutine  BODAN  calls  subroutine  BODIS  which  is  the 
routine  that  clusters  according  to  Eq.  (1).  If  the  body  is  not  analytic,  a 
table  of  x,y  ordinates  are  read  from  BODY.  These  ordinates  are  normalized  and 
then  distributed  as  a  function  of  axis  length  (chord)  or  arc  length  with  calls 
to  BODIS.  The  newly  distributed  body  points  are  interpolated  from  the  table 
of  ordinates  using  cubic  splines  (subroutine  CSPLIN)^  For  example,  ordinates 
of  y  versus  arc  length  s_  are  interpolated  to  form  y  ,  as  a  function  of  the 
distribution  arc  length  s.  Subroutine  BODY  then  returns  to  main.  At  this 
point  a  sting  and/or  cut  can  be  added  by  calls  to  STING  (i.e. ,  a  sting  as  in 
Figure  1  or  upper  cut  as  in  Figure  4)  and  CARTB  (i.e.,  forward  cut  in  Figure 
2,  lower  cut  in  Figure  4). 

Points  are  distributed  along  the  sting  and/or  cut  using  BODIS.  In  sub¬ 
routine  STING,  points  are  read  in  from  1  to  NCGRD.  NCGRD-1  points  are  added 
to  the  total  count  in  5.  Likewise  in  CARTB  a  set  of  points  along  the  cut  are 
added  to  the  previous  number.  The  final  number  of  x,y  inner  boundary  points 
is  printed  in  main. 

B.  Outer  Boundary 

MAIN  calls  subroutine  OUTER  which  forms  the  outer  boundary.  Here  the 
cubic  segments,  as  defined  by  Eq.  (4),  are  read  in  and  joined  together. 
Allowance  is  made  for  8  possible  segments.  Finally,  points  are  distributed 
along  this  boundary  as  a  function  of  arc  length.  Either  a  uniform  distribu¬ 
tion  is  used,  or  again  subroutine  BODIS  is  employed.  Interpolation  of  the 
distributed  points  is  again  obtained  by  cubic  splines,  but  the  cubic  spline 
function  is  restricted  within  an  originally  defined  segment.  Thus  in  Figure 
12b  the  cubic  spline  interpolation  is  not  carried  from  a  to  c,  but  is  carried 
a  to  b,  b  to  c,  etc.  In  this  way  the  discontinuous  corner  is  not  spline  fit. 

C.  Grid  Generation 


Finally  MAIN  calls  subroutine  ALGRD.  In  subroutine  ALGRD  the  straight 
line  ray  grid  is  formed  using  uniform  clustering  in  n.  Any  segment  of  the 
grid  between  vertical  or  horizontal  boundaries  can  then  be  regenerated  using 
subroutine  RELAX  to  obtain  an  SLOR  solution  to  Eq.  (9).  Finally,  the  grid 
lines  can  be  exponentially  reclustered  in  n 

using  Eq.  (8).  The  grid  is  then  stored  for  display  or  computational  purposes. 
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D.  Subsidiary  Subroutines 

With  storage  of  the  grid  the  program  ends.  Besides  those  subroutines 
described  above,  several  other  routines  are  called.  Subroutine  CLUST  is 

called  by  BODIS  and  this  is  the  routine  that  literally  corresponds  to  the 
distribution  function,  Eq.  (1).  Subroutines  TRIB  and  TRIP  are  routines  for 
the  solution  of  the  tridiagonal  matrix  which  must  be  inverted  in  the  succes¬ 
sive  line  overrelaxation  procedure  used  by  RELAX. 

The  TRIB  routine  is  for  conventional  tridiagonal  matrices,  the  TRIP 

routine  is  for_  periodic  tri diagonal  matrices.  Finally,  subroutine  INI PQ  is 
used  to  input  P  and  Q  of  Eq.  (9).  Use  of  these  terms  is  not  currently  recom¬ 
mended. 

Two  subroutines  are  called  from  BODAN  to  describe  the  blunt,  secant- 
ogive,  nose  projectile.  SCALC  computes  x-values  associated  with  the  nose 

cap.  The  fuse  height  is  used  to  vary  the  degree  of  bluntness.  SECANT  is  then 
called  to  provide  the  analytic  functions  used  to  compute  the  points  along  the 
remaining  body  configuration. 

Subroutine  GSPIN  is  called  only  when  a  three-dimensional  grid  is 

required.  Both  three-  and  two-dimensional  grids  are  written,  however,  the 
former  is  used  for  flowfield  computations  and  the  latter  is  used  for  plotting. 
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PROGRAM  MAIM  I  INPUT* OUTPUT*TAPE5- INPUT, TAP£6 -OUTPUT* T APE  9 , T A PE iO I 
CHECK  FOR  30  GRIO 

COMMON  JMAX*  KMAX*  JM*  KM*  N800*  J800 

COMMON  /90U0Y/  XXI100),  YYI1001*  X$  1 100 )  *  YS(100l#  SS(IOO)*  SI100) 
1  *  T I 100 I *  TSI100) 

COMMON  /COMP/  XllOOl*  YI103) 


RE AO  (5*60)  I  30*  NO*LM AX 
WRITE  (6*50)  I  30*N0*LMAX 


9ISTRI0UTE  POINTS  ALONG  INNER  BOUNOARY 
WRITE  16*701 
CALL  BOOT 
RE AO  (5*60)  NFL AG 
IF  ( NFLAG.LT •  0 )  GO  TO  AO 
RE  AO  (5*60)  NCGRO*  NCAR  T 
WRITE  (6*801  NCGRO#  NC ART 
IF  (NCGRO.GT.O)  CALL  STING  (NCGRO) 

IF  ( NCART • GT • 0 )  CALL  CARTS  ( NC ART) 

JMAX- JBOO 

JM- JMAX-1 

WRITE  (6*90)  JMAX 

IF  (NCGRO.GT.O.OR.NCART.GT.O)  SO  TO  10 
GO  TO  30 

10  WRITE  (6*100) 

00  20  J* 1  *  JMAX 

WRITE  (6*110)  J#XX( J)#TY( J) 

20  CONTINUE 
30  CONTINUE 

FORM  OUTER  BOUNOARY 
WRITE  (6*120) 

RE AO  15*60)  NSEG  S  * IOUTO 
WRITE  (6*130)  JMAX*  NSEGS 
CALL  OUTER  (NSEGS* IOUTO) 

GIRO  GENERATION 
WRITE  ( 6* l AO ) 

CALL  ALGRO  ((STOR) 


FORM  30  GRIO,  LHAX  IS  CIRCUMFERENTIAL  OIRECTION 


IF  (I30.EQ.I 

90 

STOP 

50 

FORMAT 

(  lHOi 

60 

format 

(815) 

70 

FORMAT 

( 1H1 1 

80 

FORMAT 

(  1H0  i 

90 

format 

(  1H0  i 

100 

FORMAT 

(  IHO  i 

110 

format 

(1H  , 

120 

format 

(  IHl  * 

130 

FORMAT 

(  IHO  * 

190 

format 

(  IHl  * 

ENO 

*  15* 2F 13*6) 


VALUE  OF  JMAX  *15) 

VALUES  OF  J*  X*  Y  ALONG  INNER  BOUNOARY) 


OUTER  BOUNOARY  ♦♦♦♦♦♦♦♦♦♦) 


MAIN 

MAIN 

MAIN 

main 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 


2 

3 

9 

5 

6 
7 
6 
9 

10 

11 

12 

13 

19 

15 

16 

17 

18 

19 

20 
21 
22 
23 
29 

25 

26 

27 

28 

29 

30 

31 

32 

33 
39 

35 

36 

37 

38 

39 

90 

91 

92 

93 
99 

95 

96 

97 

98 

99 

50 

51 

52 

53 
59 

55 

56 

57 

58 


GRIO  GENERATION 


MAIN 

59 

MAIN 

60 

MAIN 

61 

37 


SUBROUTINE  ALGR3  (ISTDR) 

ALGRO 

2 

COMMON  JMAX,  KMAX,  JM ,  KM,  NB03,  J300 

ALGRO 

3 

COMMON  /  BOUOY/  XX(lOO),  YY(lOO),  XStlOO),  YS(lDO),  SSIIOD),  SI  LOO ) 

ALGRO 

4 

I  ,  T  I  LOO ) •  TSIIOO) 

ALGRO 

5 

COMMON  /GRID/  X( 80,60),  Y(80,60) 

ALGRO 

6 

COMMON  /ARRAY/  A(100),  B(IDO),  CILOO),  0(100),  FI100J,  H(lOO) 

ALGRO 

7 

ALGRO 

8 

FORM  ALGEBRAIC  GRID  OR  ELLIPTIC  EQ.  GENERATED  GRIO 

ALGRO 

9 

ALGRO 

10 

READ  15,153)  KM  AX , I T E RM , IPER , NC LU 5 , IS T OR , JE LL I 

ALGRO 

11 

WRITE  (6,160)  KMAX, I  TERM , I  PER , NCL JS, ISTOR, JELLI 

ALGRO 

12 

KM-KM AX- 1 

ALGRO 

13 

READ  (5,170)  OS, OMEGA 

ALGRO 

14 

WRITE  (6,130)  OS, OMEGA 

ALGRO 

15 

ALGRO 

16 

STRAIGHT  RAY  GRIO  USED  IF  ITERM  .LE.  0 

ALGRO 

17 

OTHERWISE  US  EO  AS  INITIAL  GUESS 

ALGRO 

18 

00  40  J-  1 ,  J MAX 

ALGRO 

19 

xo-xxc  J) 

ALGRO 

20 

X1«XS( J) 

ALGRO 

21 

YO- Y  Y ( J) 

ALGRO 

22 

Yl-YS ( J) 

ALGRO 

23 

R-SQRTC ( X1-X0)**2M Y l-YO )**Z ) 

ALGRO 

24 

EPS-0* 

ALGRO 

25 

IF  ( ITERM. LT.O, )  GO  TO  10 

ALGRO 

26 

00«R/(KMAX-1) 

ALGRO 

27 

GO  TO  20 

ALGRO 

28 

10  EPS-EPS IL( R ,0* , OS, KMAX, 0.0000 2, 20, 4) 

ALGRO 

29 

OD-OS 

ALGRO 

30 

20  X(J,1)-X0 

ALGRO 

31 

Y  (  J  ,  1 )  -  Y  0 

ALGRO 

32 

TR-O. 

ALGRO 

33 

00  30  K-2 , <M AX 

ALGRO 

34 

TR-TR*DD*( l.*EP$)4* (K-2) 

ALGRO 

35 

TT-TR/R 

ALGRO 

36 

X(  J,KI  «XOMX1-XO)*TT 

ALGRO 

37 

Y(J,KI-YO*(  Yl-YOIMT 

ALGRO 

38 

30  CONTINUE 

ALGRO 

39 

40  CONTINUE 

ALGRO 

40 

ALGRO 

41 

IF  (  ITERM. _T. 01  GO  TO  70 

ALGRO 

42 

ALGRO 

43 

ELLIPTIC  P.D.E.  GRIO  GENRAT I ON  SCHEME 

ALGRO 

44 

IF  (JELLI. GE.l)  GO  TO  50 

ALGRO 

45 

CALL  RELAX  ( I  TERM , I  PER , 1 , JHA X , 3ME GA ) 

ALGRO 

46 

GO  TO  70 

ALGRO 

47 

50  00  60  LL-1, JELLI 

ALGRO 

48 

READ  (5,150)  J I , JF 

ALGRO 

49 

WRITE  (6, IDO)  J I , JF 

ALGRO 

50 

60  CALL  RELAX  ( I  TERM , I  PER , J I , JF , 3M EG  A ) 

ALGRO 

5i 

70  CONTINUE 

ALGRO 

52 

ALGRO 

53 

CLUSTERING  OPTION 

ALGRO 

54 

IF  (NCLUS.-T.O)  GO  TO  130 

ALGRO 

55 

00  120  J - 1 , J  MAX 

ALGRO 

56 

T( 11-0. 

ALGRO 

57 

00  80  K-2, KMAX 

ALGRO 

58 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


83 

EPS -EPS  IL<  T  (  XMAX  )  *  Da  *DS*KNAX,D.D0D02  *20*  J  I 
5(11-0. 

N--1 

DD  90  X-2,<MAX 
N-N^  1 

90  SIX  I -S(X-t I ^DS* ( 1. *EPSI PPN 
DO  100  X  ■! «  XMAX 
XXtXI -X( J*<l 
t DD  YY(XI-Y(J,<) 

CALL  CSPLIN  (S*XS*T*XX*A.8.C*D*F*H*2*XM*1*XNAX) 

CALL  CSPLIN  (S.Y5.TtYY,A*B,C*D*F*H*2,XM*l*<MAXI 
DD  110  X-2.XM 
X( J.KI-XSIK) 

113  Y(J,XI-YS(XI 
120  CONTINUE 

130  CONTINUE 

SUPPRESS  THE  PRINTOUT 
X  -  KMAX/2 
OD  20  J-l.JMAX 

WRITE (6*603  I  X(J*2I*Y(J*2I*X(J*3»*Y(J*3»*X(J*A>*Y(J*AI* 
1  X( J«KI  ,Y( J,X| fX( J.XMI .Y(J.XM) 

603  FDRMATdH  «1DF12.5) 

20  CONTINUE 


OPTIONAL  STORE  OF  DATA 

IF  (ISTOR.LE.O)  CO  TO  1*0 

WRITE  (91  ( (X(J*KI. J-t*JNAXl  .X-l.XNAXI *4 (VI J*X)  .J-l.JMAX) *<-l* 

1  XMAX ) 

1*0  CONTINUE 
RETURN 


150  FORMAT  (8151 

160  FORMAT  (1H0.*  SU8.  ALCRO  PRINTOUT.*..  XN AX *  I  TER N* I  PER *NCLUS * ISTOR • 
1  JELLI *  */ *8  I  51 
170  FORMAT  (8F10.0I 

180  FORMAT  (IHO.llH  OS.  OMEGA  «2F13.5) 

190  F  DR  MAT  ( 1H0 *  26HELL IPT IC  CRIO  FORMED  FROM  *I3*5H  TO  *131 
END 


ALCRO 

59 

ALCRO 

60 

ALCRO 

61 

ALCRO 

62 

ALCRO 

63 

ALCRO 

6* 

ALCRO 

65 

ALCRO 

66 

ALCRO 

67 

ALCRO 

68 

ALCRO 

69 

ALCRO 

70 

algrd 

71 

ALCRD 

72 

ALCRO 

73 

ALCRD 

7* 

ALCRD 

75 

ALCRO 

76 

ALCRD 

77 

ALCRO 

78 

ALCRD 

79 

ALCRO 

80 

ALCRD 

81 

ALCRD 

82 

ALCRD 

83 

ALCRD 

8* 

ALCRO 

85 

ALCRO 

86 

ALCRO 

87 

ALCRD 

88 

ALCRO 

89 

ALCRO 

90 

ALCRO 

91 

ALCRO 

92 

ALCRD 

93 

ALCRO 

9* 

ALCRO 

95 

ALCRO 

96 

ALCRO 

97 

ALCRO 

98 

ALCRO 

99 

ALCRO 

100 

3  0 


SUBROUTINE  B  00  AN  USEGS1 

COMMON  /CALC/  XO,  X  ATF  ,  BS,  ICT,  FLAG.  POYOX,  FXF,  RFN 
COMMON  / l  NP  S  /  XI,  X2,  X  3 »  XA,  RAO,  DYDX,  CHORO,  FUSE,  AS,  RADS 
COMMON  JMAX,  KM A  X ,  JM,  KM,  NBOO,  JBOO 

COMMON  /  BOUOY  /  XX(IQO),  YYIlOO),  XSUOO),  YSIIOO),  SSIIOO),  SHOO! 
1  ,  TI1001,  TS I 1001 

RE  AO  IS.IOD)  TAU,FLAG 
WRITE  <6, 110)  TAU,FLAG 
JWR IT- JBOO 
ICT-0 

I F  I FLAG • GE • 1 • I  CALL  SCALC 
CALL  BOO  IS  I  ISEGS, JWR  IT) 

00  10  J* l , J  BOO 
10  XXI J) -SI J) 

ANALYTIC  BODY  SHAPE 

IF  (FLAG. GE. 0.1  GO  TO  30 
00  20  J- 1 , JBOO 
XXI Jl-SI JI/SI JBOO) 

PARABOLIC  ARC 

20  YYI J)»2. ATAU^XX I J I ♦ I l.-XXI J) 1 
GO  TO  90 
30  CONTINUE 

IF  IFLAG.LE.O.I  GO  TO  AO 
CALL  SECANT 
GO  TO  90 

PROJECTILE  WITH  TANGENT  OG I VE , CYL I NOER , BOATTA I L ♦  CIRCULAR  CAP 
AO  RE  AO  IS, 100)  XI, X2, X3,XA, RAO, THETA, CHORO 
THETA-THETAA .017AS33 
OYOX-TANI TNETA) 

BASE -RAO* I XA-X3) • DYDX 
XE-BASEM  OYOX^SORTl  l.  ♦OYOXAA  2  I  I 
XS-XA*XE 

NONOIMENSIONAL  OPTION 

IF  ICHORO.SE .0. I  CHQRD-X  S-Xl 
CHORO- A8SI CHOROI 

WRITE  lb, 120)  XI  ,X2, X3.XA, RAO, THETA, CHORO 

RCH- l. /CHORO 

X1-X1ARCH 

X2-X2ARCH 

X3-X3ARCH 

XA-XAARCH 

XS-XSARCH 

BASE-BASEARCH 

XE-XE*RCH 

RAO-R AOARCH 

00  BO  J-l.JBOO 

IF  I XXI J) .GE.X2)  GO  TO  SO 

X0G-X2-XI 

XBAR-I XXI J ) — X l 1 /XDG 
VLAM-XOG/RAO 
VS3-VLAMAA2 
RBAR-VS0*.2S 

RADl-l.-VS0All.-XBAR)AA2/|RBARA*2) 

YYI J  t  * « 1.-2.ARBARAI 1,-SORT  I  RAO  I III  ARAD 
GO  TO  BO 

so  IF  I  XXI J I . GE .X3 I  GO  TO  60 
YYI JI-RAD 
GO  TO  BO 

60  IF  I  XXI J ) . GE . XA )  GO  TO  70 
YYI JI-RAD*! XXI J1-X3 I •DYDX 
GO  TO  80 

70  RS-l BASE A*2 I  At l. ♦0Y0XAA2 I 

XBAR-lXXl J)-XA) 

RA01-RS-IXBAR-BASEA0Y0X)AA2 

YYI J) -SORTIR AOl) 

80  CONTINUE 
END  OF  PROJECTILE 
90  CONTINUE 
RETURN 

100  FORMAT  IBFIQ.O) 

HO  FORMAT  IIH0,BHTAU,FLAG,2FIA.S)  runon  /  7F1A  5) 

120  FORMAT  (1H0.29H  XI  ,X2  ,X3.  XA,  KO.  TH6T  A. CHORD  ./.7FM.5I 

EM0  40 


BOOAN 

2 

BOOAN 

3 

BOOAN 

A 

BOOAN 

S 

BOOAN 

6 

BOOAN 

7 

BOOAN 

8 

BOOAN 

9 

BOOAN 

10 

BOOAN 

11 

BOOAN 

12 

BOOAN 

13 

BOOAN 

IA 

BOOAN 

IS 

BOOAN 

16 

BOOAN 

17 

900AN 

ia 

BOOAN 

19 

BOOAN 

20 

BOOAN 

21 

BOOAN 

22 

BOOAN 

23 

BOOAN 

2A 

BOOAN 

2  S 

BOOAN 

26 

BOOAN 

27 

BOOAN 

28 

BOOAN 

29 

BOOAN 

30 

BOOAN 

31 

BOOAN 

32 

BOOAN 

33 

BOOAN 

3  A 

BOOAN 

3S 

BOOAN 

36 

BOOAN 

37 

BOOAN 

38 

BOOAN 

39 

BOOAN 

AO 

BOOAN 

A  l 

BOOAN 

A2 

BOOAN 

A) 

BOOAN 

AA 

BOOAN 

AS 

BOOAN 

A  6 

BOOAN 

A  7 

BOOAN 

A  9 

BOOAN 

A9 

BOOAN 

SO 

BOOAN 

SI 

BOOAN 

S  2 

BOOAN 

S3 

BOOAN 

SA 

BOOAN 

SS 

BOOAN 

S6 

BOOAN 

S  7 

BOOAN 

sa 

BOOAN 

S  9 

BOOAN 

60 

BOOAN 

61 

BOOAN 

62 

BOOAN 

63 

BOOAN 

6  A 

BOOAN 

6S 

BOOAN 

66 

BOOAN 

67 

BOOAN 

68 

BOOAN 

69 

BOOAN 

70 

BOOAN 

71 

BOOAN 

72 

BOOAN 

73 

BOOAN 

7  A 

BOOAN 

7S 

BOOAN 

76 

BOOAN 

77 

t  • 

o  o  O 


SUBROUTINE  BOOIS  1 1 SEGSt JWRI T I 

COMMON  /CALC/  XO,  XATF,  BS*  ICTt  FLAG*  FOYOX*  FXF,  RFN 
COMMON  JMAXt  UMAX «  JMt  KM*  NBODt  JBOO 

COMMON  /BOUOY/  XXUOOlt  YYCIOOI,  XStlOOl*  YSIIOOI*  SSClOOIt 
i  •  T ( 100 ) •  TSCIOOI 

S  0 1  SIR IBUT I  ON  ON  BODY 

00  10  I  *lt ISEGS 

REAO  ( S*  20 1  JI. JFfXI .XF.OXItO XF 
WRITE  (6,301  JltJFtXlfXFfOXIfOXF 
CALL  CLUST  ( Jit jFtXltXFfOXIfOXF,SI 
10  CONTINUE 

WRITE  ( 6, AO )  (S( JltJ*l,JNRITI 
RETURN 

20  FORMAT  ( 2I*,6F10«0) 

30  FORMAT  ( 1H0 • 21H  J 1 1 JF • X I • XF, OX  I • OXF  #2I5tAF12.51 
AO  FORMAT  ( 1H  *10F11.5) 

ENO 


BOOIS 

2 

BOOIS 

3 

BOOIS 

A 

BOOIS 

5 

BOOIS 

6 

BOOIS 

7 

BOOIS 

B 

BOOIS 

9 

BOOIS 

10 

BOOIS 

11 

BOOIS 

12 

BOOIS 

13 

BOOIS 

1A 

BOOIS 

IS 

BOOIS 

16 

BOOIS 

17 

BOOIS 

IB 

BOOIS 

19 

BOOIS 

20 

BOOIS 

21 

41 


SUBROUTINE  BOOT 

BOOY 

2 

COMMON  JMAX,  KM AX «  JM,  KM*  NBOO*  JBOO 

BOOY 

3 

common  /bouoy/  xx(ioo)*  rrciooi*  xsaoo>,  rsciooi*  ssuoo),  s(ioo) 

BOOY 

A 

l  *  T(100),  TS ( 1 00) 

BOOY 

5 

COMMON  /ARRAY/  A  <  1 00 )  *  B(100>,  C(iOO),  0(100),  F(100),  H(IOO) 

BOOY 

6 

COMMON  /COM p/  X( 100 ) ,  Y(100l 

BOOY 

7 

BOOY 

a 

READ  (5*110)  NBOO, JBOO, IXORStlSEGS 

BOOY 

9 

WRITE  (6*120)  NBOO, JBOO, IXORS, ISEGS 

BOOY 

10 

IF  NBOO  IS  NEGATIVE,  ANALYTIC  SHAPE  IS  USEO 

BOOY 

n 

IF  (NBOO. GE. 01  GO  TO  10 

BOOY 

12 

NBOO— NBOO 

BOOY 

13 

CALL  BOOAN  (  ISEGS) 

BOOY 

1A 

GO  TO  90 

BOOY 

15 

BOOY 

16 

10 

CONTINUE 

BOOY 

17 

WRITE  (6*130) 

BOOY 

IB 

READ  (5,150)  CHORO 

BOOY 

19 

WRITE  ( 6 , 1 AO )  CHORO 

BOOY 

20 

RE  AO  CAROS  IN  CLOCKWISE  NOSE  TO  TAIL 

BOOY 

21 

00  20  J- 1  *  NBOO 

BOOY 

22 

READ  (5,150)  X(J)*Y(J) 

BOOY 

23 

WRITE  (6,160)  J,X(J),Y(J) 

BOOY 

2  A 

X(J)-X( J)/CHORO 

BOOY 

25 

Y( J)-Y( J )/CHORO 

BOOY 

26 

20 

CONTINUE 

BOOY 

27 

COMPUTE  NORMALIZED  ARC  LENGTH  TO  USE  AS  A  MONOTONE  PARAMETER 

BOOY 

26 

55( l)-0. 

BOOY 

29 

00  30  J" 2  *  NBOO 

BOOY 

30 

30 

5S( J)-55(  J-i)«-SQRT( ( X ( J ) —X ( J— 1 ) ) **2*( Y( Jl-Y ( J-l ) )**2) 

BOOY 

31 

00  AO  J-2*NB00 

BOOY 

32 

AO 

SS ( J ) -SS ( J ) /SS(NBOO) 

BOOY 

33 

BOOY 

3A 

COMPUTE  A  NORMALIZED  CLUSTERED  PARAMETRIC  FUNCTION  FOR  OISTRIBUTUON 

BOOY 

35 

OF  BOOT  POINTS 

BOOY 

36 

J WR I T  * JBOO 

BOOY 

37 

CALL  BOO  I S  ( ISEGS, JWRIT) 

BOOY 

38 

SAV-S(l) 

BOOY 

39 

5( L)-0. 

BOOY 

AO 

00  50  J-2, JBOO 

BOOY 

Ai 

SOEL-ABS (S( J ) -SAY) 

BOOY 

A2 

5AV-5( J) 

BOOY 

A3 

50 

5( J)-5( J-i)>SDEL 

BOOY 

AA 

00  60  J* 1 , J  BOO 

BOOY 

A  5 

60 

5( J)-5( J )/5*  JBOO) 

BOOY 

A6 

OPTION  WHETHER  TO  SET  X  EQUAL  TO  S  DISTRIBUTION 

BOOY 

A7 

IF  ( IXORS.LT. 0)  GO  TO  70 

BOOY 

A6 

CALL  CSPLIN  (S*XX*SS,X*A,B,C*0,F*H,1* JBOO, l, NBOO) 

BOOY 

A9 

CALL  CSPLIN  ( S,YY,SS,Y*A*B,C*0,F ,H,1, JBOO, 1, NBOO) 

BOOY 

50 

GO  TO  90 

BOOY 

51 

70 

CONTINUE 

BOOY 

52 

00  BO  J-l, JBOO 

BOOY 

53 

90 

XX ( J)-X( 1)>5( J)*(X(NBOO)-X(L> ) 

BOOY 

5  A 

CALL  CSPLIN  (XX,YY,X,Y,A,B,C,0,F,H*l,J9O0,l,NBO0) 

BOOY 

55 

BOOY 

56 

90 

CONTINUE 

BOOY 

57 

WRITE  (6*1701 

BOOY 

58 

00  100  J-l*  JBOO 

BOOY 

59 

WRITE  (6,160)  J,XX( J ) *YY( J) 

BOOY 

60 

100 

CONTINUE 

BOOY 

61 

RETURN 

BOOY 

62 

BOOY 

63 

110 

FORMAT  ( A  I 5  ) 

BOOY 

6  A 

120 

FORMAT  ( IHO,22HNB03,JBOO,IXORS, ISEGS  ,AI5) 

BOOY 

65 

130 

FORMAT  (1H0»38H  X,Y  INPUT  OEFINING  BOOY...  SUB.  900Y  1 

BOOY 

66 

1  AO 

FORMAT  (1H0,32H  NORMALIZING  CHORD  LENGTH  INPUT  ,F13.5J 

BOOY 

67 

150 

FORMAT  (2F10.0) 

BOOY 

66 

160 

FORMAT  (1H  * 16H  J,X,Y,  ON  BOOY  ,I5,2F1A.5) 

BOOY 

69 

170 

FORMAT  ( 1H1 , 32H  J  ANO  BOOY  DISTRIBUTED  X  ANO  Y  ) 

BOOY 

70 

FNn 

-  - 

42 


SUBROUTINE  GSPlN  ( I 30  *  NO  *  I STOR * LN AX ) 

GSPIN 

2 

COMMON  J MAX *  KMAX*  JM,  KM,  NB0O*  J8D0 

GSPIN 

3 

COMMON  /GRIO/  X  ( 80*60)*  Y(30*60) 

GSPIN 

4 

COMMON  /GR I 030 /  X3 ( 420  *60 ) *  Y3(420*60)*  23(420*60) 

GSPIN 

5 

LEVEL  2«X3*Y3*23 

GSPIN 

6 

PI-4.4ARCTANU.) 

GSPIN 

7 

OS • P  I  / ( LMAX-3 ) 

GSPIN 

8 

Ni-KMAX*LMA  X 

GSPIN 

9 

00  10  K« 1  *  KMAX 

GSPIN 

10 

OT  — 2.40S 

GSPIN 

11 

00  10  L-  1  *  -  MAX 

GSPIN 

12 

OT*OT*OS 

GSPIN 

13 

KL-(K-l) *N3*L 

GSPIN 

14 

00  10  J-1*JMAX 

GSPIN 

15 

Y3(KL»J)-Y( J*K)*SIV(DT) 

GSPIN 

16 

23(KL*J)-Y( J*K)4C0S(0T) 

GSPIN 

17 

X  3 ( KL  *  J ) -X(  J,K) 

GSPIN 

18 

10  CONTINUE 

GSPIN 

19 

IF  (ISTOR..E.O)  GO  TO  40 

GSPIN 

20 

REKINO  9 

GSPIN 

21 

WRITE  (9)  (  <  X3(  KL* J ) *KL-1 *N1 )* J-l* JMAX) ,((Y3IKL*J)*KL-1*N1),J-1* 

GSPIN 

22 

1  JMAX), <  <  23(KL*J)*KL-1*N1)* J-l* JMAX) 

GSPIN 

23 

GSPIN 

24 

REWRITE  30  DATA  FOR  20  PLOTTING 

GSPIN 

25 

GSPIN 

26 

00  30  J-1*JMAX 

GSPIN 

27 

00  20  KL-2*  Nl *LMAX 

GSPIN 

26 

N« 1 ♦ 1 KL-2 ) / NO 

GSPIN  , 

29 

XI  J*N)-X3(KL* J) 

GSPIN 

30 

Y ( J«N ) -2  3 1 KL  «  J ) 

GSPIN 

31 

20  CONTINUE 

GSPIN 

32 

30  CONTINUE 

GSPIN 

33 

WRITE  (10)  ( ( XI J*N) * J-i* JMAX) *N-1,<MAX) , 1 IYI J*N) * J-l* JMAX) *N-1* 

GSPIN 

34 

1  KMAX) 

GSPIN 

35 

GSPIN 

36 

40  CONTINUE 

GSPIN 

37 

WRITE  (6*90) 

GSPIN 

36 

00  60  L-2*LMAX*3 

GSPIN 

39 

00  60  K-i*KMAX*4 

GSPIN 

40 

KL- ( K-l ) *NQ*L 

GSPIN 

41 

WRITE  (6*70)  L*K*L*KU 

GSPIN 

42 

KL2-K*N0*L 

GSPIN 

43 

00  50  J- 1  *  J  MAX  *2 

GSPIN 

44 

50  WRITE  (6*80)  J * X3 ( KL * J  )  * V3 (KL * J ) *  231 KL * J ) « J* X3 ( KL2 * J ) * Y3 1 KL2 t J ) * 

GSPIN 

45 

1  23 ( KL2  *  J ) 

GSPIN 

46 

60  CONTINUE 

GSPIN 

47 

RETURN 

GSPIN 

46 

GSPIN 

49 

70  FORMAT  ( 1  HO  *  5M  J  *2HK-,I2*3H  L- *  1 2  * 9X * 1HX * 1DX * 1HY * 12X * 1H2 *  1 3X * 

GSPIN 

50 

1  3HJ  *2HK-*I2*3H  L-*I2*5X*1HX*10X*1HY*11X*1H2) 

GSPIN 

51 

80  FORMAT  ( 1H  *  I 3  *  1 4X * F 10. 5  * 3X* F 1 0 . 5 * 3X , F10. 5* 9X *  I 2  * 14X* F10. 5* 2 X *F10 . 

GSPIN 

52 

15*2X*F10.5) 

GSPIN 

53 

90  FORMAT  (1H1) 

GSPIN 

54 

ENO 

GSPIN 

55 
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SUBROUTINE  OUTER  ( MSECS* IOJTQ) 

OUTER 

2 

COMMON  J M  AX  *  KMA  X  ,  JM,  KM,  N803*  JBOO 

OUTER 

3 

COMMON  / BOUOY/  XX(IOO),  YY(IOO),  XS(IOO),  YS(IOO),  SS(IOO),  S(IOO) 

OUTER 

A 

1  ,  T( 100) «  TS(IOO) 

OUTER 

5 

COMMON  /COMP/  XllOOh  Y  (  100) 

OUTER 

6 

COMMON  /ARRAY/  A(130),  8(100),  C(100),  0(100),  F(100),  H(IOO) 

OUTER 

7 

2 

OUTER 

B 

C 

THIS  PROGRAM  FORMS  AN  OUTER  GRIO  BOUNOARY  USING  CONTIGUOUS  CUBIC 

OUTER 

9 

c 

SEGMENTS-  NUMBER  OF  SEGMENTS  IS  NSEGS.  POINT  ANO  SLOPE  ARE  INPUT  AT 

OUTER 

10 

c 

THE  ENOS  OF  A  SEGMENT.  SLOPE  IS  AN  ANGLE  IN  OEGRESS.  PARAMETRIC 

OUTER 

11 

c 

CuaiCS  USEO  TO  PERMIT  ANY  SLOPE  (THETA  *  90,  -90,  ETC).  INITIAL 

OUTER 

12 

c 

LOGIC  DETERMINES  CUBIC  COEFFICIENTS  OF  EACH  SEGMENT.  REMAINING 

OUTER 

13 

c 

LOGIC  DISTRIBUTES  POINTS  ALONG  OUTER  BOUNOARY  USING  ARC  LENGTH 

OUTER 

1A 

w 

AS  DISTRIBUTION  FUNCTION.  THUS  TWO  PARAMETRIC  VARIABLES  ARE  USEO. 

OUTER 

15 

r 

FINDING  X.Y  SO  CUBIC  SEGMENTS  CAN  BE  DISJOINT  IN  SLOPE  IS  MESSY 

OUTER 

16 

w 

4  SINGLE  SPLINE  INTERPOLATION  CANNOT  BE  USEO  OVER  THE  COMBINED 

OUTER 

17 

c 

SEGMENTS  BECAUSE  OF  POSSIBLE  SLOPE  DISCONTINUITY. 

OUTER 

18 

z 

i 

OUTER 

19 

DIMENSION  J A ( B ) ,  JB ( B ) 

OUTER 

20 

DIMENSION  CAO(B),  CA1 ( 6 ) ,  CA2(8),  CB0(8>,  CBl ( B ) ,  CB2(B),  CARC(8) 

OUTER 

21 

00  110  N-l.NSEGS 

OUTER 

22 

z 

POINTS  ANO  SLOPES,  -90  .LE..  THETA  .LE.  90  OEGREES  USEO 

OUTER 

23 

READ  ( 5 , 2  AO )  XO, Y0,X1, Y1 ,THO, TH1 

OUTER 

2A 

WRITE  (6,250)  X0,Y0,X1, Y1,TH0, TH1 

OUTER 

25 

RTH0-0.017A53292PTH0 

OUTER 

26 

RTH1-0.017A53292PTH1 

OUTER 

27 

XI-X1-X0 

OUTER 

2  B 

ETA-Y1-Y0 

OUTER 

29 

z 

OUTER 

30 

c 

SET  MFLAG,  LOGIC  CHIEFLY  USEO  TO  AVOIO  INFINITE  OY/OX,  USES  OX/DY-O 

OUTER 

31 

T  A*S I N (R  THO ) 

OUTER 

32 

T8-COS(RTHO ) 

OUTER 

33 

TC  *S I N( RTH1 ) 

OUTER 

3  A 

TO-COS(RTHl) 

OUTER 

35 

IF  (ABS(TA) .GT.ABS(TB))  GO  TO  13 

OUTER 

36 

NFL A  G* 1 

OUTER 

37 

IF  (ABS(TC) .GT.ABS(TO))  HFLAG-2 

OUTER 

3B 

GO  TO  20 

OUTER 

39 

10  MFLAG  *3 

OUTER 

AO 

IF  (ABS(TC) .GT.ABS(TO) )  MFLAG-4 

OUTER 

A  1 

20  CONTINUE 

OUTER 

A2 

c 

OETERMINE  COEFFICIENTS  FOR  PARAMETRIC  CUBICS 

OUTER 

A3 

c 

SET  UP  LINEAR  COEFFS.  FIRST.  INDEFINITE  CUBIC  DEFAULTS  TO  LINEAR 

OUTER 

AA 

4  l  *  X  I 

OUTER 

45 

42*0. 

OUTER 

46 

Bl-ETA 

OUTER 

47 

82*0  • 

OUTER 

48 

c 

OUTER 

49 

GO  TO  (30,40,50,60),  MFLAG 

OUTER 

50 

30  DYOXO*  T  A/ T3 

OUTER 

51 

OYOX1-TC/TO 

OUTER 

52 

TEST-0Y0X1-0Y0X0 

OUTER 

53 

IF  (ABS(TEST).LT. 0.0005)  GO  TO  80 

OUTER 

54 

SMl-2.P(ETA-0Y0X0PXI )/TEST 

OUTER 

55 

SM0*2.*XI-SM1 

OUTER 

56 

SNO*SMO*OYOXO 

OUTER 

57 

GO  TO  70 

OUTER 

58 

44 


()Ou 


40  0Y0X0-TA/T8 
0X0Y1-T0/TC 
TE  ST  *1 .—0 Y0X0*0X0Y1 
IF  (A0S( TEST). LT. 0.00051  G3  TO  80 
SM0-2.M  XI-0X0Y1*ETA J/TEST 
SN0-SM0*0Y0X0 
GO  TO  70 
53  0X0Y0-T8/TA 
0Y0X1-TC/T0 
TEST-1. -0Y0X1*0X0Y3 
IF  (A0S(TEST).LT. 0.00051  GO  TO  80 
SN0-2.META-0Y0X1AXI  J/TEST 
SHO-SNOAOXOYO 
GO  TO  70 
60  0X0Y0-T8/TA 
0X0Y1-T0/TC 
TEST-0X0Y1-0X0Y0 

IF  (ABS(TEST).LT. 0.0005)  GO  TO  80 
SN1-2.M  XI-OXOYOAETAJ/TEST 
SN0-2.AETA-SN1 
SMO-SNOAOXOYO 

73  Al-SMO 
A2-XI-SM0 

81- SNO 

82- ETA-SNO 
80  CONTINUE 

4N8R-25 

COMPUTE  NUMERICAL  ARC  LENGTH  AS  A  PARAMETER!  EXACT  ARC  LENGTH  IS 
POSSIBLE  8UT  INVERSE  PROCESS  IS  NOT 
00  90  J- 1  *  JN 8R 
T  T  *  4N8R-1 
TT-( J-ll/TT 
TUJ-XO*TTM  A1^A2*TT) 

90  Y( 4)-Y0*TTM81*82ATT) 

S  ARC-0. 

NOTE  ...COULO  USE  SUN  OP  SQUARES  AS  PARAMETER  RATHER  THAN  ARC  LENGTH 
IN  THIS  WAY  ONE  CAN  AVOIO  SQUARE  ROOT  CALCULATION  MUST  USE  EVERY 
WHERE 

00  100  4-2 1 JN8R 

100  SARC-SARC+SORTl ( X ( J ) -X ( 4-1 I) *P2* ( Y ( 4 >-Y < 4-1 I |**2) 

OATA  FOR  EACH  CU8 I C  SEGMENT 
CAl (N)-Al 
CA2(N)-A2 
CBKNI-81 
C82CNI-82 
CAOINI-XO 
CBO(N|-YO 
CARCf  Nl -S ARC 

WRITE  (6*  260 )  XO * Ai, A2, YO, 01 ,02 ,SARC 
110  CONTINUE 


CU3ICS  OETERMINEO,  NOW  OISTRI0UTE  POINTS 


OUTER 

59 

OUTER 

60 

OUTER 

61 

OUTER 

62 

OUTER 

63 

OUTER 

6* 

OUTER 

65 

OUTER 

66 

OUTER 

67 

OUTER 

68 

OUTER 

69 

OUTER 

70 

OUTER 

71 

OUTER 

72 

OUTER 

73 

OUTER 

7* 

OUTER 

75 

OUTER 

76 

OUTER 

77 

OUTER 

70 

OUTER 

79 

OUTER 

80 

OUTER 

81 

OUTER 

82 

OUTER 

83 

OUTER 

as 

OUTER 

85 

OUTER 

86 

OUTER 

07 

OUTER 

60 

OUTER 

89 

OUTER 

90 

OUTER 

91 

OUTER 

92 

OUTER 

93 

OUTER 

9* 

OUTER 

95 

OUTER 

96 

OUTER 

97 

OUTER 

90 

OUTER 

99 

OUTER 

100 

OUTER 

101 

OUTER 

102 

OUTER 

103 

OUTER 

10  A 

OUTER 

105 

OUTER 

106 

OUTER 

107 

OUTER 

100 

OUTER 

109 

OUTER 

110 

OUTER 

111 

OUTER 

112 

OUTER 

113 

OUTER 

114 

OUTER 

115 

45 


TOTAL  OUTER  ARC  LENGTH 

OUTER 

OUTER 

116 

117 

SARC-O. 

OUTER 

1  IB 

00  120  N-1,NSEGS 

OUTER 

119 

120  SARC-SARC*CARC(N) 

OUTER 

120 

WRITE  (6*270)  SARC 

OUTER 

121 

OEFINE  A  UNIFORM  OUTER  DISTRIBUTION  ARC  LENGTH 

OUTER 

OUTER 

122 

123 

RH-l./( JMAX-1) 

OUTER 

124 

00  130  J* 1  *  J  MAX 

OUTER 

125 

130  SS(J)«( J-1)*RH 

OUTER 

126 

OPTIONAL  USE  OF  CUBIC  SEGMENTS  T3  CLUSTER 

OUTER 

OUTER 

127 

128 

if  ( i out o* le • 0)  go  to  iso 

OUTER 

129 

call  BOOIS  (IOUTO,JNAX) 

OUTER 

130 

00  1 AO  J-1*JNAX 

OUTER 

131 

l AO  SS( J)-( S ( J)-S( 1) )/(  S< JHAXI-SC 1)  ) 

OUTER 

132 

150  CONTINUE 

OUTER 

133 

WRITE  (6*280)  (SS( J)* J-i* JMAX) 

OUTER 

13  A 

NORMALIZE  OUTER  ARC  LENGTH  SEGMENTS  TO  SCALE  OF 

OISTRIBUTION 

ARC  LEN  OUTER 

135 

CA-O. 

1  OUTER 

136 

00  160  N-1,NSEGS 

OUTER 

137 

CA-CA*CAfcC(N)/SARC 

OUTER 

138 

160  CARC(N)-CA 

OUTER 

139 

WRITE  (6*290)  ( C ARC ( N ) * N- 1* NS E GS ) 

OUTER 

1  AO 

FINO  J  INDICES  LIMITS  WITHIN  A  SEGMENT 

OUTER 

1A1 

N-i 

OUTER 

142 

JA(N)-1 

OUTER 

1A3 

00  180  J-2*JHAX 

OUTER 

lA  A 

IF  (SS( J).LE.CARC(N) )  GO  TO  173 

OUTER 

145 

N-N*l 

OUTER 

146 

JA(N)-J 

OUTER 

147 

170  J  8 ( N ) * J 

OUTER 

1A8 

180  CONTINUE 

OUTER 

149 

00  190  N-l,NSEGS 

OUTER 

150 

WRITE  (6*300)  J A ( N ) *J8(N) 

OUTER 

151 

190  CONTINUE 

OUTER 

152 

FORM  PARAMETRIC  AR  R  AT  S  *  FROM  OISTRIBUTEO  PARAMETRIC  ARRAY* 

OUTER 

OUTER 

153 

154 

JSE  IT  TO  OETERMI NE  X*Y  WITHIN  A  OUTER  SEGMENT 

CURVE. 

OUTER 

155 

SPLINE  REQUIRES  ABOUT  5  POINTS  IN  AN  INTERVAL 

OUTER 

156 

S(l)-0. 

OUTER 

157 

RT-l./( JNBR-l) 

OUTER 

158 

00  220  N-UNSEGS 

OUTER 

159 

T ( 1 ) *0* 

OUTER 

160 

IF  ( N*  GT • 1 )  S( l)-CARC(N-l) 

OUTER 

161 

X(1)-CA0(N) 

OUTER 

162 

Y( i)-CBO(N) 

OUTER 

163 

00  200  J-2*JNBR 

OUTER 

164 

TT-(  J-1)'*RT 

OUTER 

165 

TU)-TT 

OUTER 

166 

X( J ) -CAO( N) ♦TTA(CAi ( N)*TT*CA2(N) ) 

OUTER 

167 

Y(  J)*CB0(N)-*TT*(CB1(N)*TT*C82(  M )  ) 

OUTER 

16B 

OS*SQRT( (X( J)-X( J-l) )4*2MY( J  )  -  Y (J-l) )**2) 

OUTER 

169 

S( J ) -S ( J-l)*OS/SARC 

OUTER 

170 

WRITE  (6*280)  T(J)*X(J)*Y(J).S(J) 

OUTER 

171 

200  CONTINUE 

OUTER 

172 

46 


210 

220 

230 

2*3 

250 

260 

270 

2  80 

290 

300 

310 


173 

17* 

175 

176 

177 

178 

179 

180 

181 

182 

183 

18* 

185 

186 

137 

188 

139 

190 

191 

192 

193 

19* 

195 

196 

mi 


Jl» JA (N) 

J2» J81NI 

CALL  CSPLIN  tSS*TS*S*T*A*8*C*0tF* 
00  210  J-J1.J2 
TT-TSl JJ 

XSt  J)-CA0tNH>TT*(CAHN|*TT*CA2tNI 

TSt J)-C80<M)*TT*CC3ltNI^TT*Ca2tNl 

WRITE  (6*280)  TSC JI*XS( 4)*YS< J) 

CONTINUE 

CONTINUE 

00  230  J* 1*  JM  AX 

WRITE  (6*310)  XS(J)*YS(J) 

CONTINUE 

RETURN 


H*  Ji* J2*l* JNBR ) 


) 

) 


FORMAT  (8F10.0) 

FORMAT  (1H0,21H  XO *  TO  *  X 1  * T 1  * THO • THi  «6F13*5) 

FORMAT  ( 1H0*23HX0*A1*A2*Y0*81*82*SARC  */*7Fl3.5) 

FORMAT  ( 1H0  «  *2  H  OUTER  80UN0AR Y  NONO I  ME  NS ! ONAL  ARC  LENGTH  *F1**6) 
FORMAT  (1H  «  8F 12 • 5 ) 

FORMAT  (1H0.6H  CARC  *6F13.5) 

FORMAT  (1H  *7H  JA,J8  *215) 

FORMAT  (1H  • 6HXS  «  YS  «2F15.6) 

P  MD 


OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 

OUTER 
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SUBROUTINE  RELAX  C I  TERM*  IPER , JI , JF ,OMEGA ) 

RELAX 

2 

COMMON  JMAXt  KMAXt  JMt  KM,  N83D,  JBOO 

RELAX 

3 

COMMON  /GRIO/  X( 90,60 I ,  YC80,50I 

RELAX 

4 

COMMON  /ARRAY/  AC130I,  8C100I,  CC100I,  OCIOOI,  FC1D0I,  HCXOOl 

RELAX 

5 

OIMENSION  IPCl42if  IRC142I 

RELAX 

6 

0 I  ME  NS  I  ON  GC100I 

RELAX 

7 

COMMON  /SOURCE/  P  C  80 , 2  1 1  OC8Dt2l,  PFAC( 21  *  0FACI2I 

RELAX 

8 

RELAX 

9 

SLOR  SOLUTION  OF  ELLIPTIC  GRIO  GENERATION  EOS. 

RELAX 

XO 

EL  XI  ANO  OEL  ETA  -  1.0 

RELAX 

XX 

RELAX 

12 

Jl-l 

RELAX 

13 

J2-JMAX 

RELAX 

14 

IF  CIPER.GT.OI  GO  TO  10 

RELAX 

X  5 

ji*jm 

RELAX 

16 

J2-JF-X 

RELAX 

X  7 

ID  CONTINUE 

RELAX 

18 

CALL  INIPO 

RELAX 

19 

RELAX 

20 

ITER-0 

RELAX 

21 

KM-KHAX-l 

RELAX 

22 

SET  PERIOOIC  INOICES 

RELAX 

23 

00  20  J- 1 , J MAX 

RELAX 

24 

IP( JI* J*l 

RELAX 

25 

20  IRC  J  1  - J— 1 

RELAX 

26 

IPC  JMAX) -1 

RELAX 

27 

IRC ll-JMAX 

RELAX 

28 

RELAX 

29 

FORM  DIFFERENCE  EXPRESSIONS  ANO  TRIOIACONALS 

RELAX 

30 

30  I TER- I TER* X 

RELAX 

3  X 

RSUM-O. 

RELAX 

32 

00  X 60  K K -2 , KM 

RELAX 

33 

K-<m*2-KK 

RELAX 

34 

KP-K 

RELAX 

35 

KR-K-1 

RELAX 

36 

CP  1-EXP  C  (X-KIPPFACCin 

RELAX 

37 

CP2-EXPC U-KMAX  I  *PF  ACC  21 I 

RELAX 

38 

cox  -exp  c  ci-k i*of  acc  in 

RELAX 

39 

C02-EXPC  CK-KMAX) P0FACC2) 1 

RELAX 

40 

00  *0  J- J 1 , J2 

RELAX 

*1 

JP* I P C  J 1 

RELAX 

42 

JR- 1 R  C JI 

RELAX 

43 

XXO-CXC JP,< 1-XC  JR,< 1 |p,5 

RELAX 

44 

XEO-CXC  J,KPI— XC  J,KR) )*•* 

RELAX 

45 

YXO-C YC  JP,K 1-YC  JR,<) I  ♦*  5 

RELAX 

46 

YEO-CYC J.KPl-YC J,KR1 1*,5 

RELAX 

67 

A0-XE0**2*YE0**2 

RELAX 

48 

8  0— 2.  PC  XXO*XEO*YXO*YEOI 

RELAX 

49 

G0-XX0**2*YX0*P2 

RELAX 

50 

XXE0-.25«C  XC  JP«KPI-XC  JP* KR l-XC JR«KP1 ♦XC J*,K*1 1 

RELAX 

5  X 

YXE0-.25M  YC  JP ,KP 1 -Y C  JP , KR l-Y ( JR , K P I ♦ Y C  J* , XR I  I 

RELAX 

52 

AC  JI-AO 

RELAX 

53 

8 ( J 1 -- AD- AD- GO-GO 

RELAX 

54 

CC JI-AO 

RELAX 

55 

FC  JI—  BOPXXEO-GOMXC  J,KPI*XC  Jf  <R  1 1 

RELAX 

56 

GC  J I  — 80PYXE0-G0PC  Y C  J,KP 1 ♦YC J,<RI  1 

RELAX 

57 

SOURCE  TERMS 

RELAX 

58 

o  o  o  f  i  n  n 


OJAC»(XXO*YEO-XE  0*YXO) 

0SQ»OJAC4*2 

COFA»OSOMCP1*P  <  Jf  1DCP2*P<  Ji  2)  > 

COF9-OSQ*! CO 1*Q< Jf LI»C02PQ(  J»2) ) 

F0BX-SIGNI.5.C0F A) 

FOBE-SIGNl .5.COFB) 

A(  J)»A(  JI-COFAM.5-FOBXI 
B( J|-B( J)-2.*ICOFA*FOBX*COFB*FOBE) 

Cl J)-Cl J)*C0FA*l.5*F0BX) 

F(  J  )«Fl  JI-COFBM  l.5»F0BE )*Xl J,<P)-t .5-FOBE)*Xl  J»KR)) 
Gt  J)»GI J)-COFB*t (.5»FOBEl*Y( J, <P l-< . 5-FOBE 1*Y<  J*KRII 

hq  continue 

IF  IIPER.GT.O)  GO  TO  L30 
SET  B.C.  AND  INVERT 

XI  UN  ANO  1  AX  PLANES  MUST  BE  X  3R  Y  CARTESIAN  PLANES 
OUTFLOW  B  C  ON  X 

TEST  WHETHER  XI  PLANE  IS  X  OR  Y  •  CONSTANT  PLANE 
IF  (ABS  ( Y (  JF»KNAX)-Y(  JF*  LU-LT. O.OOll  GO  TO  50 
:  OUTFLOW  XI-PLANE  TAKEN  AS  X-CONSTANT  PLANE 

A ( JF  )»0. 

B<  JF)  —  L. 

CUFJ-O. 

FIJFI—  XI  JF«1I 
GO  TO  60 

C  OUTFLOW  XI-PLANE  TAKEN  AS  Y«  CONSTANT 
50  A(JF)-1. 

3 I  JF)  —  L. 

Cl JFI-O. 

FUF)-(-XUF— itKDXI JF-2*K) )/3. 

60  CONTINUE 
c  INFLOW  B  C  ON  X 

IF  (ABSIYI JItKMAX)-Y(JI»l)).LT.  0.001)  GO  TO  70 
C  INFLOW  XI-PLANE  TAKEN  AS  X  -  CONSTANT  PLANE 
A( JI )-0. 

B(  JI)— L. 

CUII-O. 

FIJD-XIJhl) 

GO  TO  BO 
70  CONTINUE 

C  INFLOW  XI-PLANE  TAKEN  AS  Y  -  CONSTANT  PLANE 
A( JI )-0. 

BUI)— L. 

CUIl-1. 

FIJI )• I— X IJ I ♦ 1»K ) ♦<  U I tK ) 1 /3. 

90  CALL  TRIB  I  A « B «C «0« F « J I « JF) 

C  Y  B.C. ANO  INVERSION 

IF  I  ABS  I Y  UF  »KHAX) -Y(  JF*  1 )  ).LT. 0.001  )  GO  TO  90 
C  OUTFLOW  XI-PLANE  TAKEN  AS  X-CONSTANT  PLANE 
A I  JFI-1. 

B(  JF)—  1. 

C( JF )-0. 

G( JF)-l-Y(JF-ltK) +r I JF-2.K) )  7  3. 

GO  TO  100 
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62 

RELAX 

63 

RELAX 

64 
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67 

RELAX 

6B 

RELAX 

69 

RELAX 

70 

RELAX 

71 

RELAX 

72 

RELAX 

73 
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7* 

RELAX 

75 

RELAX 

76 

RELAX 

77 

RELAX 

7B 

RELAX 

79 

RELAX 

BO 

RELAX 

B 1 

RELAX 

62 

RELAX 

S3 

RELAX 

64 

RELAX 

B5 

RELAX 

B6 

RELAX 

B7 

RELAX 

B6 

RELAX 

89 

RELAX 

90 

RELAX 

91 

RELAX 

92 

RELAX 

93 

RELAX 

94 

RELAX 

95 

RELAX 

96 

RELAX 

97 

RELAX 

96 

RELAX 

99 

RELAX 

100 

RELAX 

101 

RELAX 

102 

RELAX 

103 

RELAX 

104 

RELAX 

105 

RELAX 

106 

RELAX 

107 

RELAX 

108 

RELAX 

109 

RELAX 

110 

RELAX 

111 

RELAX 

112 

RELAX 

113 

RELAX 

114 

RELAX 

115 
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C  OUTFLOW  XI-PLANE  TAKEN  AS  Y  •  CONSTANT  PLANE 
90  ACJFI-O. 

at  jf  >  — 1. 

C( JF )*D. 

GCJFI— YOFtll 
C  INFLOW 

100  IF  IA8SIT* J I*KNAX)-Y( Jl*l ) ).LT.0.00i)  GO  TO  HO 
C  INFLOW  XI-PLANE  TAXEN  AS  X  -  CONSTANT  PLANE 
AC JII-O. 

B( 

C( JI 1-1. 

G(JI)*-fT(JI»l»<)-r(Jl>2»K))/3. 

GO  TO  120 

C  INFLOW  XI-PLANE  TAXEN  AS  Y  *  CONSTANT  PLANE 
110  ACJD-O. 

3 ( JI  )*-L. 

Cl  JII-O. 

G(  JII—  V(  JI*1> 

120  CALL  T R I  8  ( A  *  8 *C *  0. G . J I  * JF ) 

GO  TO  140 

c 

130  CONTINUE 

C  PERIOOIC  a.c. 

CALL  TRIP  ( A*B*C*F*D*H*l * JMAX) 

CALL  TRIP  ( A*B*C*C*0*H*1* JHAX) 

140  CONTINUE 

RELAXATION  UPOATE 
00  150  J - J I  *  JF 
YC-GC  J  >- Y ( J*X) 

XC-FC  J)-X( J.K) 

X( JtK)-X( J*  X  I ♦OHEG A4XC 
YC J*K)-YC J*K»^OHEGAPYC 
150  RSUM-RSUN+ABS(XC)>ABS(rC) 

160  CONTINUE 

IF  (( ITER/101410. LT. ITER)  GO  TO  170 
WRITE  (6*180)  ITERtRSUH 
170  CONTINUE 

IF  ( I  TER .LT • I  TERN )  GO  TO  30 
RETURN 

180  FORMAT  ( 1H  *26H  ITERATION  NBR  ANO  RSUN  *I5*E12.3) 
END 


RELAX  116 

RELAX  117 

RELAX  118 

RELAX  119 

RELAX  120 

RELAX  12  L 

RELAX  122 

RELAX  123 

RELAX  124 

RELAX  125 

RELAX  126 

RELAX  127 

RELAX  128 

RELAX  129 

RELAX  130 

RELAX  131 

RELAX  132 
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RELAX  134 

RELAX  135 
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RELAX  137 

RELAX  138 

RELAX  139 

RELAX  140 

RELAX  141 

RELAX  142 

RELAX  143 

RELAX  144 

RELAX  145 

RELAX  146 

RELAX  147 

RELAX  148 

RELAX  149 

RELAX  150 

RELAX  151 

RELAX  152 

RELAX  153 

RELAX  154 

RELAX  155 

RELAX  156 

RELAX  157 

RELAX  158 

RELAX  159 

RELAX  160 
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SUBROUTINE  INIPQ 

COMMON  J M AX «  K  MAX  *  JM*  KM*  NSOO*  JSOO 
COMMON  /SOURCE/  P(80*2)*  0<80*2>*  PFAC(2)» 
DIMENSION  PC ( 2  )  *  0C(2I 


QF AC  I  2) 


READ 
RE  AO 
WRITE 
WRITE 


(3*20)  PFAC(1),QFAC(1)*PFAC(2),0FAC<2) 

(5*20)  PC(1)»3C(1)*PC(2).3C(2) 

(6*30)  PF AC(1)*0FAC( 1 ) * PF AC ( 2 ) » OF  AC  1 2 
(6*40)  PC(1)*OC(1)*PC(2)»OC(2) 


) 


10 


00  10  N*  1  *  2 
00  10  J* 1  *  J  MAX 
P(  J*N) -PC(N) 

Q( J  *N)-QC (N ) 

return 


20  FORMAT  (8F1Q.Q) 

30  FORMAT  ( 1H3  *  *  EXPONENT  COEFF 
1  • * / *4F 13*  5  ) 

40  FORMAT  ( 1H0  * 18H  PC1*0C1*PC2* 
ENO 


I  Cl E  NTS  FOR  SOURCE 
QC2  f  AF 1 3«  5 ) 


TERMS»PFAC*OFAC 


SUBROUTINE  CLUST  ( J I  * J F *  X I *XF *0 XI  * DXF * S ) 
0  I  ME  NS  I  ON  S(l) 

XFXI-XF-XI 
H  *  1 « / ( JF— J I ) 

H2*H*H 

H3*H2*H 

C-(0XF»0XI-2.*H*XFXI )/ (H— 3*PH2^2«*H3) 

8 ■ ( OX  I  — H*XF  XI-CMH3-H) )/ (H2-H) 

A-XFX I-8-C 

00  10  J  ■  J I  *  JF 
X-( J-JI)*H 

io  s (  j)-xi>xm k+x*i 

RETURN 

ENO 
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RELAX 
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RELAX 

164 

RELAX 

165 
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RELAX 

178 

RELAX 

179 

RELAX 

180 

RELAX 
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RELAX 
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RELAX 

183 

WORKS 

2 

WORKS 

3 

WORKS 

4 

WORKS 

5 

WORKS 

6 

WORKS 

7 

WORKS 

8 

WORKS 

9 

WORKS 

10 

WORKS 

11 

WORKS 

12 

WORKS 

13 

WORKS 

14 

WORKS 

15 

WORKS 

16 

WORKS 

17 
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WORKS 

16 

FUNCTION  EPSIL  I F NX  * FM IN , OF* * NPT * FPCC *  I CC * NC ALL 1 

WORKS 

19 

WORKS 

20 

THIS  SUBROUTINE  APPLIES  A  NEWTON-RAPHSON  ROOT-FINOING 

WORKS 

21 

TECHNIQUE  TO  FIND  A  VALUE  OF  EPSILON 

for  a  PARTICULAR  USE 

WORKS 

22 

OF  THE  EXPONENTIAL  STRECHING  TRANSFORMATION. 

WORKS 

23 

WORKS 

24 

FNX  IS  TOTAL  ARC  LENGTH  ALONG  COOROINATE 

WORKS 

25 

FMIN  IS  STARTING  VAL'JE  OF  ARC  LENGTH  SUCH  AS  0.0 

WORKS 

26 

OF N  IS  SPECIFIEO  INITIAL  INCREMENT  OF  ARC 

LENGTH 

WORKS 

27 

NPT  IS  NUMBER  OF  POINTS  ALONG  C03R0INATE 

WORKS 

26 

FPCC  IS  ITERATIVE  ERROR  BOUNO,  E.G.OI  0.000021 

WORKS 

29 

ICC  IS  MAXIMUM  NUM8ER  OF  ITERATIONS 

WORKS 

30 

NC  ALL  IF  NCALL«1  INITIAL  GUESS  FOR  EPS 

IS  USEO 

WORKS 

31 

IF  NC  ALL  .GT.  1*  PREVIOUS  EPS  USEO  AS  INITIAL  GUESS 

WORKS 

32 

WORKS 

33 

WORKS 

3* 

fmxl-fmk 

WORKS 

35 

fminl-fmin 

WORKS 

36 

OFNL-OFM 

WORKS 

37 

FPCCL-FPCC 

WORKS 

38 

ICCL-ICC 

WORKS 

39 

WORKS 

40 

WORKS 

41 

FNPTM2-NPT-2 

WORKS 

42 

IF  INCALL. EQ.l)  E P S- I FMKL/OF ML  1 ♦ ♦ 1 1. 0/FNPTM2 l-l • 0 

WORKS 

43 

WORKS 

44 

00  10  NIT-1, ICCL 

WORKS 

45 

EPl-EPS^i.3 

WORKS 

46 

EP1TN-EP1**FNPTM2 

WORKS 

47 

REPS-i. O/EPS 

WORKS 

48 

DFMOE-OFML*REPS 

WORKS 

49 

F-FMXL-FN I NL-OFHOE  ♦lEPiTNPGPi-i.Ol 

WORKS 

50 

IF  I  ABSIFl.LT. FPCCLI  CO  TO  20 

WORKS 

51 

0FM0E2-0FM0E*REPS 

WORKS 

52 

FPN-0FM0E2*!  1.0^EPlTN*IEPSPFNPTH2-i.0n 

WORKS 

53 

EPS-EPS^F/FPN 

WORKS 

54 

10  CONTINUE 

WORKS 

55 

WORKS 

56 

EPSIL-EPS 

WORKS 

57 

WRITE  <6,301 

WORKS 

56 

RETURN 

WORKS 

59 

WORKS 

60 

20  EPSIL-EPS 

WORKS 

61 

SUPPRESS  THE  ESPIL  PRINTOUT 

WORKS 

62 

HR  I TE I  6 , 601 1  EPSIL, F, NIT 

WORKS 

63 

WORKS 

64 

RETURN 

WORKS 

65 

WORKS 

66 

WORKS 

b7 

30  FORMAT  I/42H  EKCEEOEO  MAX.  NO.  OF  ITERATIONS  IN  EPSIL. 1 

WORKS 

66 

ENO 

WORKS 

69 
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c 

c 

c 

c 


c 

c 

c 


c 

c 


c 

c 


c 


WORKS 

70 

SUBROUTINE  CSPLIN  (XX*YY*X,Y  ,A,8,C,0,F,H,N1 

,N2, J1,J2» 

WORKS 

71 

DIMENSION  XX  ( 1 )  t  YYCl),  XCl),  YCl),  A(l),  BCD,  C(l>*  Dll)*  F(l>, 

WORKS 

72 

1  NC  1) 

WORKS 

73 

WORKS 

IS 

CUBIC  SPLINE  INTERPOLATION 

WORKS 

75 

X,Y  ARRAYS  ARE  TO  BE  INTERPOLATED 

WORKS 

76 

YY  ARE  FOUND  INTERPOLATES  CORRESPONDING  TD  XX 

WORKS 

77 

J 1 tt  J2  ,  ARE  INOICE  LIMITS  ON  X,Y 

WORKS 

78 

Nl,  N2  ARE  INDICE  LIUTS  ON  XX  (  ALSO  YY) 

WORKS 

79 

DIMENSION  OF  ARRAYS  CARRIED  IN  FRDM  OUTSIDE,  XCJ)  MUST  BE  MONOTONE 

WORKS 

80 

ORMULA  FROM  NUMERICAL  METHODS  BY  DA HLQU I S T , 8 JORCK , A NOER SON 

WORKS 

81 

JLS  FEB.  77 

WORKS 

82 

WORKS 

83 

ROUNDING  ERROR  PROTECTION  CAUTION  MAY  MASK  ERRDR  IN  LOGIC 

WORKS 

8* 

IF  (XXCND.LT. XC  Jl)  >  XXCND-XCJl) 

WORKS 

95 

IF  (XXCN2I.GT.XC J2) )  XXCN2)-XCJ2) 

WORKS 

86 

WORKS 

87 

WORKS 

88 

FIRST  FIND  DERIVATIVE  LIKE  TERMS  THAT  ARE 

COEFFICIENTS 

WORKS 

89 

J  A  ■  J  1  ♦  1 

WORKS 

90 

J8- J2-1 

WORKS 

91 

00  LO  J-JA,J2 

WORKS 

92 

H(  J)-X( J)-X(  J-l) 

WORKS 

93 

10  0(J)"(Y(J)-Y(J—1))/H(J) 

WORKS 

9* 

00  20  J  *  J  A  , J8 

WORKS 

95 

A C  J )*H( J ♦ 1) 

WORKS 

96 

8C  J ) *2, *  C HC J  DH( J ♦ 1 >  ) 

WORKS 

97 

20  F(  J )-3.+(H(  J)*0(  J*n*H(4n»A0UM 

WORKS 

98 

B( Jl)-2. 

WORKS 

99 

H  (  J  1  >  -  I  . 

WORKS 

100 

F  (  JD-3.P0C  JA) 

WORKS 

101 

A(J2)-1. 

WORKS 

102 

8C J21-2. 

WORKS 

103 

F ( J  2 ) *  3 • *0( J2) 

WORKS 

IQS 

CALL  TRIB  (A,8,H,C,F,J1, J2 ) 

WORKS 

105 

WORKS 

106 

INTERPOLATION  ,  XCJ)  ARRAY  MUST  BE  MONOTONE 

WORKS 

107 

J  *  J 1 

WORKS 

108 

I  * ji ♦! 

WORKS 

109 

00  80  N* N l , N2 

WORKS 

110 

30  IF  CXC J).LE.XX(N) .ANO.XC I)#SE.XX(N) )  GO  TO 

70 

WORKS 

111 

IF  (X(I)-XX(N))  <»0,  **0 , 50 

WORKS 

112 

*0  J-J*l 

WORKS 

113 

I  -  I  ♦  1 

WORKS 

IIS 

IF  (1.GT.J2)  GO  TO  60 

WORKS 

115 

GO  TO  30 

WORKS 

116 

50  J*J-1 

WORKS 

117 

I-I-i 

WORKS 

118 

IF  CJ.LT.J1)  GO  TO  60 

WORKS 

119 

GO  TO  30 

WORKS 

120 

60  WRITE  (6,90) 

WORKS 

121 

STOP 

WORKS 

122 

70  T-CXXCN) -XC J) )/H( I) 

WORKS 

123 

TT-l.-T 

WORKS 

12* 

YY(N)«T*Y(  I  )*TT*YC  JDHI  l  )*T*TTM  IFC  J  )-DU  1  )*TT-I  F(  I  )-0(  I  )  )*T) 

WORKS 

125 

BO  CONTINUE 

WORKS 

126 

RETURN 

WORKS 

127 

WORKS 

128 

90  FORMAT  C1HD,20H  ERROR  IN  CSPLIN  ) 

WORKS 

129 

ENO 

WORKS 

130 
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SUBROUTINE  TR!B  ( A, B • C, X ,F ,NL, NU ) 

0  ME  NS  I  ON  A<2>,  B(2),  C<2),  X<2>,  FC2) 
XML)*C(NL)/B*NL) 

F ( NL  )*FI NL) /B(NL) 

NLPL-NL+L 
00  10  J*NLP1iNU 
Z-1./(B( JJ-AI JMXU-1U 
XIJJ-CIJ»*Z 

10  F I J )•( FC J)-A( JJ*F ( J-l) MZ 
NUPNL*NU*NL 
00  20  J 1  *NLP 1 ,NU 

j-nupnl-ji 

20  F(  J)*F(  J  J-X  (JI*F(M) 

RETURN 


SUBROUTINE  TRIP  <  A,  B  ,  C  ,  F  ,0,  S  ,  J 1  ,  J2  ) 

DIMENSION  A  (  3  )  ,  9(3),  C(3),  F ( 3  )  ,  Q ( 3  )  ,  S<3> 

J  A*  J 1*1 
FN«Ft  J2) 

C  FORWARD  ELIMINATION  SWEEP 
0(  Jl)— C(  Jl)/B<  Jl) 

F  (  J 1 ) *F ( Jl) /B ( J1 ) 
sun-Aun/Bui) 

00  10  J-JA,J2 

p - 1 - / c  a i j ) ♦ a i j)*o<  j-i) ) 

0 ( J ) — C ( J  ) *  P 

F  (  J  )  » (  F  ( J  )  -  A  ( J)*F(J-1))*P 
S(  Jl—  A(  J  )  *St  J-l )*P 
10  CONTINUE 
C  BACKWARD  PASS 
J  J  *  J  1  ♦  J2 
0( J  2 ) *0. 

S ( J2I-1. 

00  20  I • J A, J2 
J-J J-I 

S(J)«S(J)OlJI*S(J*l) 

20  0  (  J ) *F ( J 1 I JI*0(J*1I 

F( J2)*(FN-C( J2)*0(J1)-AI J  2 ) • 0 ( J2-1) )  / < C < J 2 ) *S < Jl ) *A < J2 ) *S  IJ2- 1 ) 

i  ♦  auzn 

C  8ACKWAR0  ELIMINATION  PASS 
00  30  I-JA,J2 
J*J J-I 

30  F(J)*F(J2)*S(JM0( J) 

RETURN 

6N0 
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WORKS 

132 

WORKS 

133 

WORKS 

13* 

WORKS 

135 

WORKS 

136 

WORKS 

137 

WORKS 

138 

WORKS 

139 

WORKS 

1*0 

WORKS 

1*1 

WORKS 

1*2 

WORKS 

1*3 

WORKS 

1** 

WORKS 

1*5 

WORKS 

1*6 

WORKS 

1*7 

WORKS 

1*8 

WORKS 

1*9 

WORKS 

150 

WORKS 

151 

WORKS 

152 

WORKS 

153 

WORKS 

15* 

WORKS 

155 

WORKS 

156 

WORKS 

157 

WORKS 

158 

WORKS 

159 

WORKS 

160 

WORKS 

161 

WORKS 

162 

WORKS 

163 

WORKS 

16* 

WORKS 

165 

WORKS 

166 

WORKS 

167 

WORKS 

1 6  B 

WORKS 

169 

WORKS 

170 

WORKS 

171 

WORKS 

172 

WORKS 

173 

WORKS 

17* 

WORKS 

175 

WORKS 

176 

WORKS 

177 
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SUBROUTINE  SCALC 

COMMON  /CALC/  XO  ,  XATF,  BS,  ICT,  FLAG,  POYDX,  F  XF,  RFN 
COMMON  / 1 NP  S /  XI,  X2,  X3,  X*,  RAD,  DYOX,  CHORO,  FUSE,  AS,  RAOS 

THIS  SUBRTN  CALCULATES  VALUES  NEEDED  IN  THE 
EQUATIONS  SOLVING  FOR  Y  VALUES  IN  SBRTN  BODAN 

R  E AO  (5,101  X1,X2, X3,X*, RAD, THETA, CHDRO 

READ  (5,101  RADS, FUSE, AS 

THETA-THETA*. 017*533 

DYDX-TAN( THETAI 

CHORO-ABS (CHORDI 

FUSE-FUSE/2 

TO  FIND  THE  Y  VALUE  (BSI  OF  CIRCLE  USED  IN 
THE  SECANT  OGIVE  CA LC JL A T I D NS 

XBAR - AS-X2 

YS-RADS**2-XBAR*#2 

YSS-SORT(YS) 

BS-RAD-YSS 


TO  FIND  THE  X  VALUE  ( XATF I  AT  THE  FUSE 

YR-ABS(BSI*FUSE 
YRR-RA0S**2-YR**2 
XSS-SOR  T ( YRR I 
XATF-AS-XSS 


TO  FINO  THE  SLOPE 

POYOX«(XATF-AS  J/(FUSE-BS I 
POYDX— POYOX 

TO  FIND  THE  X  VALUE  (XOl  AT  THE  NOSECAP 

XSO-SORT( 1. ♦POYDXI 
RFN-FUSEPXSO 
FXF2-RFN**2-FUSE**2 
F  XF-SGRT ( FXF2 ) 

XS-RFN-FXF 

XO-XATF-XS 

RRITE  (6,201  X1,X2,X3,X*, RAD, THETA, CHORD 
HRITE  (6,301  RADS, FUSE, AS, BS,XO, XATF 
RETURN 


10  FORMAT  (8F10.0I 

20  FORMAT  ( 1H0,27HX1,X2,X3, X*, RAD, THETA, CHORD, /,«F1*. 5) 
30  FORMAT  ( 1H0 , 23HR AOS , FUSE , AS, BS, XO, XATF , /6F1*. 5  I 
END 


SCALC 

2 

SCALC 

3 

SCALC 

* 

SCALC 

5 

SCALC 

6 

SCALC 

7 

SCALC 

8 

SCALC 

9 

SCALC 

10 

SCALC 

11 

SCALC 

12 

SCALC 

13 

SCALC 

1* 

SCALC 

15 

SCALC 

16 

SCALC 

17 

SCALC 

IS 

SCALC 

19 

SCALC 

20 

SCALC 

21 

SCALC 

22 

SCALC 

23 

SCALC 

2* 

SCALC 

25 

SCALC 

26 

SCALC 

27 

SCALC 

28 

SCALC 

29 

SCALC 

30 

SCALC 

31 

SCALC 

32 

SCALC 

33 

SCALC 

3* 

SCALC 

35 

SCALC 

36 

SCALC 

37 

SCALC 

38 

SCALC 

39 

SCALC 

*0 

SCALC 

*1 

SCALC 

*2 

SCALC 

*3 

SCALC 

** 

SCALC 

*5 

SCALC 

*6 

SCALC 

*7 

SCALC 

*8 

SCALC 

*9 

SCALC 

50 

SCALC 

51 
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SUBROUTINE  SECANT  SECANT  2 

COMMON  /CALC/  XO,  X  ATF  ,  BS,  ICT,  FLAG,  POYOX,  F  XF  ,  RFN  SECANT  3 

COMMON  /INFS/  XI,  X2,  X  3 ,  X4,  RAO,  OTOX,  CHORD,  FUSE,  AS,  RAOS  SECANT  A 

COMMON  J MAX ,  X HA  X ,  JM,  KM,  NB30,  JBOO  SECANT  5 

COMMON  /SOUOY/  XX(IOO),  YY(IOO),  XS(IOO),  YSllOOl,  SS(IOO),  S(IOO)  SECANT  6 

1  ,  T(IOO),  TSI100I  SECANT  7 

C  SECANT  8 


C  PROJECTILE  WITH  NOSE CA P  ,  SE CANT  OGIVE 
C 

RCH-1. /CHORD 
X 1 ■ X 1 ♦RCH 
X2-X2PRCH 
X  3»  X  3  ARCH 
X  4*  X4*RCH 
RAO-RAO^RCH 
F  USE  »FUS  E  *R  CH 
RAOS-RAOS*RCH 
X  ATF  *X  AT  F *RCH 
AS  * AS  ARCH 
9  S  »B  S ♦RCH 
00  AO  J>I, JBOO 
IF  f  XX ( J ) «  GE • X ATF )  GO  TO  10 
C 

C  COMPUTE  Y  VALUES  FOR  NOSECAP 
C 

XB  AR  *X  X  (  J  )  -(  X  ATF  ♦FXF  > 
RA0I»RFN**Z-X8AR**2 
YYUI-SORTIRAOI  I 
GO  TO  AO 

10  IF  (XXI Jt,GE*X2t  GO  TO  20 
C 

:  COMPUTE  Y  VALUES  FOR  OGIVE 

r 

XBAR-XXI JI-AS 
RA3IS-RA0S**2-XBAR**2 
YYI J 1 -BS ♦SQRTIRAOIS I 
GO  TO  AO 

20  IF  (  XX (  J ) • GE • X3 )  GO  TO  30 
C 

C  COMPUTE  Y  VALUES  FOR  CYLINDER 
C 

YYCJ)*RAO 
GO  TO  AO 

c 

C  COMPUTE  Y  VALUES  FOR  3 OATT  A  I L 

r 

30  YY( J)«RA0+(XX(J)-X3)*0Y0X 
AO  CONTINUE 

WRITE  (6,50)  POYOX,RFN,FXF 
RETURN 

C 

50  FORMAT  ( LH0,13HP0Y0X,RFN,FXF,5X,3F 
ENO 


CYLINDER, BOATTAIL  SECANT  9 

SECANT  10 

SECANT  11 

SECANT  12 

SECANT  13 

SECANT  IA 

SECANT  15 

SECANT  16 

SECANT  17 

SECANT  18 

SECANT  19 

SECANT  20 

SECANT  21 

SECANT  22 

SECANT  23 

SECANT  2 A 

SECANT  25 

SECANT  26 

SECANT  27 

SECANT  28 

SECANT  29 

SECANT  30 

SECANT  31 

SECANT  32 

SECANT  33 

SECANT  3 A 

SECANT  35 

SECANT  36 

SECANT  37 

SECANT  38 

SECANT  39 

SECANT  AO 

SECANT  Al 

SECANT  A2 

SECANT  A3 

SECANT  aa 

SECANT  A  5 

SECANT  A6 

SECANT  47 

SECANT  48 

SECANT  49 

SECANT  50 

SECANT  51 

SECANT  52 

A. 5)  SECANT  53 

SECANT  54 


t 
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SUBROUTINE  STING  MCGROI 

STING 

2 

COMMON  JM  AX  ,  KM  AX  ,  JM,  KM,  N0OD, 

JBOO 

STING 

3 

common  /eouor/  xxuoo),  yy(ioo). 

XS(IOO),  YS(IOO),  SSI100), 

SHOO) 

STING 

4 

1  ,  T I  100 )  ,  TS(IOO) 

STING 

5 

STING 

6 

STING  ASSUMEO  ALIGNED  WITH  X  COORDINATE 

STING 

7 

STING  INPUT  AS  A  SEPARATE  0OOY,  J* 

l, NCGRO. ..BOOIS  OATA  STARTS 

AT  1 

STING 

B 

REAO  (5,20)  ISEGS 

STING 

9 

WRITE  (6,30)  ISEGS 

STING 

10 

REAO  (5,40)  XMIN,YMIN, XMAX, YMAX 

sting 

11 

WRITE  (6,50)  XMIN,YMIN,XMAX,YMAX 

STING 

12 

CALL  BOOIS  ( ISEGS, NCGRO) 

STING 

13 

00  LO  JJ-2, NCGRO 

STING 

14 

J»JBOO*J J-l 

STING 

15 

S( J J)-<S( JJ )-S( 1) ) /(S(NCGRO) -S  ( 1 

)) 

STING 

16 

XX ( J)-XMIN*S( J J )♦( XMAX-XMIN) 

STING 

17 

10  YY( J)*YMIN+S(JJ)4(YMAX-YMIN) 

STING 

IB 

JBOO*  JBOO*NCGRO- 1 

STING 

19 

WRITE  (6,60)  ( XX ( J) , J*i, JBOO ) 

STING 

20 

RETURN 

STING 

21 

STING 

22 

20  FORMAT  < 15) 

STING 

23 

30  FORMAT  (1H0,35H  STING  PART  OF  PROGRAM . ISEGS  -,I5) 

STING 

24 

40  FORMAT  (BF10.0) 

STING 

25 

50  FORMAT  (IH0,21H  XMI N, YM I N, XMAX, YMAX  ,4F13.5) 

STING 

26 

60  FORMAT  < 1H  ,8F12.5) 

STING 

27 

ENO 

STING 

28 

STING 

29 

SUBROUTINE  CARTB  (MCART) 

sting 

30 

COMMON  J  MAX ,  KMAX ,  JM  *  KM,  N0OD, 

JBOO 

STING 

31 

COMMON  /BOJOY/  XX(LOO),  YY(LOO), 

XS(IOO),  YS(IOO),  SS(IOO), 

S(100) 

STING 

32 

1  ,  T(LOO),  TS(IOO) 

STING 

33 

STING 

34 

REAO  (5,30)  ISEGS 

STING 

35 

WRITE  (6,40)  ISEGS 

STING 

36 

REAO  (5,50)  XMIN,YMIN, XMAX, YMAX 

STING 

37 

WRITE  (6,60)  XMIN,YMIN, XMAX, YMAX 

STING 

38 

SHIFT  POINTS 

STING 

39 

00  10  JS-i,JBOO 

STING 

40 

J*JBOO*NCART— JS 

STING 

41 

J  J* JBOO* i- JS 

STING 

42 

XX( J)-XX( JJ ) 

STING 

43 

10  YY(J)-YY(JJ) 

STING 

44 

J0OO »JBOO* N CART-1 

STING 

45 

CALL  BOOIS  ( ISEGS, NCART) 

STING 

46 

NCM-NCART-i 

STING 

47 

00  20  J-1,NCH 

STING 

48 

S( J)-(S( J)-S( i))/(S(NCART)-S(l) ) 

STING 

49 

XX( J)*XMIN*S( J)4( XMAX-XMIN) 

STING 

50 

20  YY( J)-YMIN*S( J ) ♦ ( YM A X-YM I N ) 

STING 

51 

RETURN 

STING 

52 

STING 

53 

30  FORMAT  (15) 

STING 

54 

40  FORMAT  ( 1H0, 33H  CARTB  PART  OF  PROGRAM . ISEGS* ,15) 

STING 

55 

50  FORMAT  (BFIO.O) 

STING 

56 

60  FORMAT  ( LHO , 21 H  XMI N ,YMI N, XMAX, YMAX  ,4F13.5) 

STING 

57 

ENO 

STING 

-5B 
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APPENDIX  B 


DEFINITION  OF  INPUT  VALUES 


Because  the  computer  code  has  so  many  options,  input  statements  are 
scattered  throughout  the  code.  Listed  below  are  the  input  parameters  to  the 
program  along  with  an  explanation  of  each  parameter.  An  input  card  is  indi¬ 
cated  below  by  numbering  and  underlining.  All  input  formats  are  either  15  or 
F10.0  and  are  so  indicated  on  the  right  hand  side  of  the  list  of  input  para¬ 
meters.  Variables  names  follow  conventional  FORTRAN  conventions  and  all 
integer  names  begin  with  I,  J,  K,  L,  N,  or  M.  Special  instructions  as  to 
whether  or  not  a  card  is  read  are  indicated  with  $  symbols. 
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INNER  BOUNDARY 


1 3D,  ND,  LMAX 

*315* 

I3D  = 

1  generates  3-dimensional 
generated. 

data,  otherwise,  2D  data  only  is 

ND  = 

same  as  LMAX. 

LMAX  = 

the  number  of  planes  in  the 

circumferential  direction. 

NBOD,  JBOD, 

IXORS,  ISEGS 

*415* 

NBOD  =  number  of  ordinates  (i.e.,  x,y  data  points)  used  to  define 
body.  If  NBOD  <  0,  an  analytic  body  shape  is  used  in  sub¬ 
routine  BODAN. 

JBOD  =  number  of  points  user  will  distribute  on  body  surface. 

IXORS...  use  x  or  s  (arc  length)  as  monotone  clustering  parameter,  x 
used  if  and  only  if  (iff)  IXORS  0. 

ISEGS  =  number  of  contiguous  clustering  segments  along  body  surface. 

Each  segment  requires  end  points  and  spacing  specification  as 
read  in  below. 

$  Read  cards  3,  4,  5  iff  NBOD  .GE.0$ 

3.  CHORD _ *F10.0* 

All  x,y  data  is  normalized  (i.e.,  divided  by)  CHORD.  CHORD  may  be  set 

to  1 .  . 

4.  X(J),  Y(J) _ *2F10.0* 

x,y  ordinates  that  define  body,  J  =  1,  NBOD  data  cards  are  read  in. 

If  x,y  are  correctly  normalized,  set  CHORD  =  1. 

5.  JI,  JF,  XI,  XF,  DXI ,  DXF _ *2I5,4F10.0* 

Data  that  defines  the  cubic  stretching  function,  see  Eqs.  (1)  and  (2) 

and  Fig.  (11) . 

There  are  ISEGS  such  cards  read-in. 
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In  the  notation  of  the  text 


JI  =  J0 
JF  =  jf 
XI  =  x0 
XF  =  xf 
DXI  =  Ax 

o 

DXF  =  Vxf 

Note  DXI  and  DXF  may  both  be  positive  or  both  be  negative  as  x  is 
increasing  or  decreasing 

$  Read  cards  6,  7,  8  iff  NBOD  .LT.O  $ 


6. 

TAU,  FLAG 

*2F10.0* 

TAU  =  parabolic  arc  thickness  ratio 

FLAG  =  0  tangent-ogive  cylinder,  boattail 

projectile  read  in 

1  secant- ogive  cylinder,  boattail 
read  in 

projectile  with  nosecap 

iff  FLAG  =  1,  card  7  goes  after  card  8b. 

7. 

JI,  JF,  XI,  XF,  DXI,  DXF 

*215 ,4F10.0* 

see  card  5 

$  iff  FLAG  ,GE.0$ 

8a. 

XI,  X2 ,  X3 ,  X4 ,  RAD,  THETA,  CHORD 

*6F10.0* 

XI  =  value  of  x  at  nose 

X2  =  value  of  x  at  ogive-cylinder  juncture 

X3  =  value  of  x  at  cyl inder-boattai 1  juncture 

X4  =  value  of  x  at  boattail  base 

RAD  =  radius  of  cylinder 

THETA  =  angle  of  degrees  that  boattail  makes  with  cylinder 
(THETA  is  negative) 

CHORD  =  if  CHORD  .GE.O,  body  length  normalized  to  one. 


Note: 


a  spherical  cap  is  added  to  boattail  so  the  body  length  is 
not  X4-X1. 
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8b.  RADS,  FUSE,  AS 


*3F10.0* 


$  iff  FLAG  . EQ. 1$ 

RADS  =  radius  of  secant 
FUSE  =  fuse  height  (at  nosecap) 

AS  =  value  of  x  at  secant-origin 

9.  NFLAG _ *15* 

Option  to  exit  program  after  body  clustering  data  is  printed  out. 
NFLAG  .LT.O,  STOP 

10.  NCGRD  ,  NCART _ *215* 

Parameters  that  allow  addition  of  sting/rear  cut  and  a  front  cut. 

NCGRD  .GT.O,  NCGRD  points  added  for  rear  cut  or  sting 

NCART  .GT.O,  NCART  points  added  for  front  cut  (or  lower  cut  of  C- 
grid) 

$  Read  cards  11,  12,  13  iff  NCGRD  .GT.0$ 

11.  ISEGS _ *15* 

ISEGS  of  sting 

12.  XMIN ,  YMIN  ,  XMAX,  YMAX _ *4F10.0* 

XMIN  =  initial  x  value  of  sting 

YMIN  =  initial  y  value  of  sting 

XMAX  =  final  x  value  of  sting 

YMAX  =  YMIN 

13.  JI,  JF,  XI,  XF,  DXI ,  DXF _ *2I5,4F10.0* 

see  card  5 


$Read  cards  14,  15,  16  iff  NCART  .GT.0$ 

14.  ISEGS  _ *15* 


15. 

XMIN,  YMIN,  XMAX,  YMAX 

*4F10.0* 

16. 

JI,  JF,  XI,  XF,  DXI,  DXF 

*215 ,4F10.0* 

Front  cut  data  like  sting  data 
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17. 


NSEGS,  IOUTD 


*215* 


NSEGS  =  number  of  contiguous  cubic  segments  which  are  used  to  form  an 
outer  boundary 

IOUTD  =  number  of  clustering  segments  along  outer  boundary  (i.e., 
previous  ISEGS) 

18.  XO,  YO,  XI,  Yl,  THO,  TH1 _ *6F10.0* 

x,y,6  end  point  values  used  to  define  cubic  segment  according  to  Eq. 
(4).  There  are  NSEGS  such  cards  read-in.  Here  0  implies  initial 
point,  1  implies  final  end  point.  The  angle  0  is  in  degrees,  and  is 
defined  in  the  usual  way.  See  discussion  of  Fig.  12  for  examples 
of  0. 

$  Read  cards  19  iff  IOUTD  .GT.0$: 

19.  JI,  JF,  XI,  XF,  DXI ,  DXF _ *215 ,4F10.0* 


See  card  5,  there  are  IOUTD  such  data  cards  read  in. 

Arc  length  clustering  used  and,  as  the  total  arc  length  is  not  known 
on  the  first  run  of  the  program  it  is  output.  Use  of  normalized  arc 
length  allows  user  to  cluster  without  true  value  of  arc  length. 


GRID  GENERATION 

20.  KMAX ,  ITERM,  IPER,  NCLUS,  ISTOR,  JELLI _ *615* 

KMAX  =  number  of  points  in  n-di recti  on 

ITERM  =  number  of  iterations  used  to  relax  Eq.  (9) 

If  ITERM  .LT.O,  straight  ray  grid  is  generated 

IPER....  set  IPER  .GT.O  if  periodic  grid  generated 

NCLUS...  NCLUS  .LT.O  means  grid  is  not  reclustered  using  Eq.  (8) 

ISTOR...  store  grid  on  computer  disc  storage  if  ISTOR  .GT.O 

JELLI...  If  JELLI  .GE.l,  Limits  JI  and  JF  are  set  on  the  elliptic  grid 
domain. 


21.  PS,  OMEGA _ *2F10.0* 

DS  =  Asq  (i.e.,  As  in  n  direction  at  n  =  0  boundary).  See  Eq. 

(8).  Note  asq  used  along  entire  n  =  0  boundary. 

OMEGA  =  relaxation  factor  for  SLOR  in  Subroutine  RELAX.  Typical  safe 
value  is  1.55.  0  <  OMEGA  <  2.0. 
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*215* 


22.  JI,  JF 

JELLI  such  cards  read  in.  Limits  of  S  min  to  S  max  over  which  an 
elliptic  solver  is  used.  JI  and  JF  must  correspond  to  vertical  or 
horizontal  rays. 

23.  BLANK  CARD _ _ _ 

24.  BLANK  CARD  _ _ _ 

"P  and  "Q  input  data,  not  recommended 
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APPENDIX  C 


SAMPLE  INPUT  AND  OUTPUT 


The  following  computer  output  illustrates  the  output  from  a  sample  grid 
generation.  The  tubular  projectile  illustrated  in  Figure  15  (the  outer  bound¬ 
ary  is  identical  to  that  shown  in  Figure  13a)  was  used  as  a  sample  case.  This 
particular  case  is  the  most  difficult  to  set  up  as  it  requires  the  largest 
number  of  special  instructions.  Input  values  are  printed  after  they  are  read 
in,  so  the  output  also  supplies  the  user  with  an  example  of  the  data  input 
cards. 
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APPENDIX  D 


PLOT  PROGRAM  LISTING 


A  listing  of  the  computer  code  used  to  generate  grid  plots  is  presented 
in  this  appendix.  The  plot  program  is  written  in  standard  FORTRAN  IV  and  uses 
the  Tektronix  Plot  10  software  package.  All  plots  were  produced  on  the 
Tektronix  4010-1  display  terminal  which  was  connected  to  the  BRL  Cyber  173/76. 
The  program  is  an  interactive  plotting  routine  which  prompts  the  user  for  all 
requested  information. 
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PROGRAM  IGP 


76/76 


OPT-L  ROUNO  • 


FTN  4.8*518 


L 


5 


LO 


15 


20 


25 


30 


35 


40 


45 


PROGRAM  IGP  ( INPUT, OUTPUT* TAPES -INPUT, TAPE  6-OUTPUT, TAPE  8 *PL3,KB I 
COMMON  /GRIOC/  GXMIn,  GXMAX,  GYMIN,  GYMAX,  Xli00,50),  Y(100,50) 
COMMON  /HOR/  TI TLE ( 5 ) 


c#####cr IP-PROGRAM  to  PLOT  THE  COMPUTATIONAL  GRIO  ABOUT  AN 
C  AEROOYNAHIC  BOOT  OR  AIRFOIL  SECTION ♦♦♦♦♦♦ 

C 

C*****INPUT  SECTION*#*#* 

C*****CONTROL  PARAMETER  INPUT 

C**NAOO--l  SUPPRESSES  LISTING  OF  PLOT  OATA,NADD-L  ALLOWS  IT 
C**JMAX,<MAX-NUM8ER  OF  J,K  POINTS 
C**HEAO— ALPHANUMERIC  INFORMATION  DESCRIBING  THE 
C  CONFIGURATION  BEING  PLOTTED 

CALL  CONNEC  (5LTAPE5) 

CALL  CONNEC  (5LTAPE6) 

CALL  TERM  ( 1 v 1 02  H ) 

CALL  SET9UF  (3) 

CALL  ANMOOE 

C  R  E  AO ( 5  9 1 0 ) HE  AO 

/• 

WRITE  (6913) 

RE  AO  (59*)  180 
180*1  BO/ 10 
CALL  INITT  (  I  BO ) 

CALL  BINITT 


C 


WRITE  (6920) 

R  E AO  (59*)  JM AX , KMAX 

WRITE  (6930) 

RE  AO  ( 5 ♦ ♦  )  GXMIN9GXMAX,GYMIN9GYMAX 
WRITE  (6,40) 

RE  AO  (5950)  TITLE 
CALL  NEWPAG 


C*  **R  E  AO  IN  X  ANO  Y  VALUES 

RE  AO  (8)  (  (X(  J,K),  J-I9  JMAX)9<-l9KMAX)  9(  (  Y  (  J  9*  )  9  J*l,  JM  AX  )  fK*  L9K.MA 
C#****  PLOT  THE  GRIO  ♦♦♦♦♦ 

CALL  GRDPLT  ( JMAX 9 < MA X9 I  BO) 

C*  *♦ TERM  I  NATE  PLOTTING 

CALL  FINITT  (0,700) 

STOP 

C 


10 

FORMAT 

(  l  9h 

WHAT 

IS 

20 

FORMAT 

(20H 

WHAT 

ARE 

30 

FORMAT 

(  30H 

WHA  T 

ARE 

*0 

FORMAT 

(  28H 

ENTER 

n 

50 

FORMAT 

END 

(5A10 

) 

3AU0  RATE?) 

JMAX,KMAX? ) 
XMIN,XMAX,YM(N,YMAX?) 
TLE  -  UP  TO  50  CHAR) 


7C 


SUBROUTINE  grdplt 


7 5/?5  DPT-l  RDUNO-^-P/ 


FTN  4.8*51B 


1 


5 


10 


15 


20 


25 


30 


35 


<•0 


SUBROUTINE  GRDPLT  I JMAX ♦ KM AX * I  30 ) 

COMMON  /GRIOC/  GXMIN*  GXNAX*  GYMIN*  GYMAX,  XI100*50)* 
COMMON  /HOR/  TITLEI5) 

DIMENSION  G  X  (  1 2  8  )  *  GYI128) 

C 

C***READJUST  PLOT  LIMITS  IN  ORDER  TO  AVOID  STRETCHED  PLOTS 
ICOUNT-0 
10  XMAX-GXMAX 
XMIN-GXM IN 
YMAX-GYMAX 
YMIN-GYM I N 
XDIF-XMAX-XMlN 
YOIF-YMAX-YMIN 
IF  IXDIF.LT.YOIFI  GO  TD  20 
XDIFH-XDIFP0.5 
Y MI D-< YMAX*YHINJ  *0*5 
YMX-YMID^XDIFH 
YMN-YMIO-XDIFH 
YMAX-YMX 
YMIN-YMN 
GO  TO  30 

20  YOI FH- YD  I F*0  •  5 

xmid-ixmax*xmin)*o*5 

XMX-XM ID*YOIFH 
XMN-XMIO-YOIFH 
X  M  A  X  -  X  MX 
XMI N-XMN 
30  CONTINUE 
C 

C  PLOT  THE  LINES 

IF  I  ICDUNT.GT.O)  GO  TO  40 
AXMlN-XMtN 
AXM  A  X *  XM  A  X 
AYMIN-YMIN 
AYMAX-YMAX 
40  CONTINUE 

c 

CALL  BINITT 
CALL  NPTS  UMAX! 

CALL  ANNODE 

WRITE  16*120)  TITLE 

CALL  XFRM  12) 

CALL  YFRM  12) 


45 


50 


55 


C 


C 


CALL 

OLI 

MX 

I  XMI  N 

*  XM AX  ) 

CALL 

OLI 

MY 

I  YMIN 

*  YMA  X ) 

CALL 

SLI 

MX 

1150* 

800) 

CALL 

SLI 

MY 

I 50,700) 

OD  70 

K- 

1* 

<M  AX 

DO  50 

J- 

1* 

J  M  A  X 

GXI  J) 

-XI 

J* 

K) 

GYC  J) 

-YC 

J  * 

<) 

IF  (K.GT.l)  GO  TO  50 
CALL  CHECK  <GX,GY) 


YI 100*50) 


79 


SUBROUTINE  GROPLT 


76/7*) 


DPT-1  ROUNO-*-*/ 


ftn  *.a*5ie 


CALL 

OSPLAY  C  GX  «  GY) 

60 

CALL 

CPLOT  ( GX*  GY) 

60 

70 

CONTINUE 

CALL 

NPTS  (KNAX) 

00  90 

J-l* JMAX 

00  00 

K-1*KMAX 

GXCK  J 

*X( J*K) 

65 

00 

GY(K  I 

-Y( J*K) 

CALL 

CPLOT  CGX,GY) 

90 

CONTI 

NUE 

CALL 

0ELL 

CALL 

TSENO 

70 

CALL 

TINPUT  (II) 

IF  C  I 

I.EQ.93)  GO  TO  110 

75 


00 


05 


90 


j 

95 


too 


CALL 

NEWPAG 

CALL 

TSENO 

00  100  J«1*JMAX 

GX(  J) 

«X(  J*l) 

GY (  J  ) 

«Y( J,l) 

CALL 

3INITT 

CALL 

OL I  NX  ( AXM (N, AXNAX ) 

CALL 

OLIMY  ( AYMIN* AYHAX) 

CALL 

SLIMX  (150*000) 

CALL 

SLIMY  (50*700) 

CALL 

NPTS  (JMAX) 

CALL 

XFRM  (2) 

CALL 

YFRM  (2) 

CALL 

CHECK  (GX,GY) 

CALL 

OSPLAY  ( GX*  GY ) 

CALL 

TSENO 

CALL 

BELL 

CALL 

ANMDOE 

WRITE 

(6,130) 

call 

TSENO 

Call 

RECOVR 

CALL 

VCURSR  (ICH,XX,YY) 

GXMIN 

l-XX 

GYMIN 

l-YY 

CALL 

ANMOOE 

WRITE 

(6*1^0  ) 

CALL 

RECOVR 

CALL 

TSENO 

CALL 

ANHOOE 

CALL 

VCURSR  (  I CH  *  XX , Y Y 1 

GXW  AX 

-XX 

105 


GYWAX-YY 


CALL  TSENO 
I  COUNT -ICOJNT*! 
CALL  NEWPAG 
GO  TO  10 

UO  110  CONTINUE 

CALL  NEWPAG 
RETURN 


120 

FORMAT 

(  5  A  10 

) 

15 

130 

1*0 

FORMAT 

FORMAT 

ENO 

(  30H 
(  IHO, 

POSITION  CURSOR  FOR  XMIn,YMIN) 

30H  POSITION  CURSOR  FOR  XMAX,YNAX) 
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this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


